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(2)  SCIENTIFIC  PROGRESS  DURING  REPORTING  PERIOD 


During  this  period  we  successfully  executed  all  elements  of  our  original  statement  of 
work.  The  work  at  Boston  University  has  been  carried  out  in  collaboration  with  Dr.  T. 
Bahder  at  AMRDEC  in  Huntsville,  Alabama. 

Our  research  program  during  this  period  has  concentrated  on  investigating  novel 
scientific  and  technological  solutions  for  ultra-sensitive  quantum-optical  registration  and 
sensing  of  phase  shifts  induced  by  external  fields.  The  use  of  several  types  of 
nonclassical  states  of  light  and  development  of  novel  ultra-sensitive  interferometric  phase 
measurement  configurations  has  been  at  the  center  of  our  research  effort. 

The  major  technical  findings  have  been  reported  in  several  leading  refereed  journal 
papers  and  presented  at  several  technical  meetings. 

The  following  results  have  been  obtained  during  the  period  of  performance: 

1)  We  performed  theoretical  investigation  of  the  limitations  that  a  single-mode 
assumption  carries  when  using  Fock-state  and  NOON-state  light  in  practical  quantum 
interferometry.  Optical  dispersion  is  a  significant  contributor  to  a  wave  packet 
modulation  in  existing  materials  demonstrating  good  phase  shifting  capability  in  response 
to  the  action  of  an  external  field.  This  degrades  the  quality  of  observed  interference 
pattern  and  the  accuracy  of  phase  evaluation. 

We  considered  Mach-Zehnder  and  Hong-Ou-Mandel  interferometers  with  nonclassical 
states  of  light  as  input,  and  studied  the  effect  that  dispersion  inside  the  interferometer  has 
on  the  sensitivity  of  phase  measurements.  We  examined  in  detail  a  number  of  different 
one-  and  two-photon  input  states,  including  Fock,  dual  Fock,  NOON  states,  and  photon 
pairs  from  parametric  down  conversion.  We  computed  the  probabilities  of  the  possible 
measurement  outcomes  as  a  function  of  the  phase  shift  introduced  between  the  arms  of 
the  interferometer,  and  then  computed  the  Shannon  mutual  information  between  the 
phase  shift  value  and  the  measurement  outcomes.  This  provides  a  means  of  quantitatively 
comparing  the  utility  of  various  input  states  and  various  interferometric  arrangements  for 
determining  phase  in  the  presence  of  dispersion.  In  addition,  we  considered  a  simplified 
model  of  parametric  down  conversion  for  which  probabilities  can  be  explicitly  computed 
analytically,  and  which  serves  as  a  limiting  case  of  the  more  realistic  down  conversion 
model. 

The  results  obtained  in  this  work  have  been  summarized  in  a  Physical  Review  A 
publication: 

•  D.  S.  Simon,  A.  V.  Sergienko,  and  T.  B.  Bahder,  ‘’Dispersion  and  Fidelity  in 
Quantum  Interferometry”,  Physical  Review  A  v.  78,  053829  (2008). 


2)  We  continued  our  quest  to  discover  alternative  quantum  techniques  for  ultra-precise 
phase-shift  measurement.  A  key  component  of  this  effort  involved  studying  methods  for 
controlling  the  dispersive  effects  that  can  degrade  the  sensor’s  effectiveness.  We 
proposed  a  novel  design  for  a  phase-sensitive  quantum  interferometer  that  demonstrates 
an  odd-order  dispersion  cancellation  effect  based  on  frequency-anticorrelated  entangled 
photons.  We  used  this  special  interferometer  to  demonstrate  observation  of  both  odd- 
order  and  even-order  dispersion  cancellation  effects  in  one  single  experiment  for  the  first 
time.  In  addition  to  demonstrating  a  peculiar  quantum  effect,  this  result  carries  strong 
promise  for  significantly  enhancing  resolution  when  evaluating  phase  shifts  in  dispersive 
optical  materials. 

These  results  provided  the  basis  for  a  research  paper  recently  published  in  Physical 
Review  Letters: 

•  Olga  Minaeva,  Cristian  Bonato,  Bahaa  E.  A.  Saleh,  David  S.  Simon,  and 
Alexander  V.  Sergienko,  “Odd-  and  Even-Order  Dispersion  Cancellation  in 
Quantum  Interferometry”,  Physical  Review  Letters,  v.  102,  100504  (2009). 

Work  is  currently  underway  for  the  experimental  investigation  of  this  idea  and  its 
application  to  quantum  phase  sensing. 

3)  During  the  course  of  our  detailed  studies  of  quantum  interference  and  phase  sensitivity 
we  were  able  to  demonstrate  a  link  between  frequency  dispersion  and  its  spatial 
counterpart,  aberration.  We  developed  a  theoretical  framework  for  the  effect  and  by 
means  of  quantum  interferometry  we  produced  the  first  experimental  demonstration  of 
even-order  aberration  cancellation.  The  effect  is  a  spatial  counterpart  of  spectral  group 
velocity  dispersion  cancellation,  which  is  associated  with  spectral  entanglement.  It  is 
manifested  in  temporal  interferometry  by  virtue  of  the  multi-parameter  spatial- spectral 
entanglement.  Spatially-entangled  photons  generated  by  spontaneous  parametric  down 
conversion  were  subjected  to  spatial  aberration  introduced  by  a  deformable  adaptive 
mirror  that  modulates  the  wavefront.  We  showed  that  only  odd-order  spatial  aberrations 
affect  the  quality  of  quantum  interference. 

Two  manuscripts  summarizing  these  findings  have  been  published  in  Physical  Review 
Letters  and  Physical  Review  A: 

•  Cristian  Bonato,  David  S.  Simon,  Paolo  Villoresi,  and  Alexander  V.  Sergienko, 
"Multiparameter  Entangled-State  Engineering  Using  Adaptive  Optics",  Physical 
Review  A  v.  79,  062304  (2009) 

•  Cristian  Bonato,  Alexander  V.  Sergienko,  Bahaa  E.  A.  Saleh,  Stefano  Bonora,  and 
Paolo  Villoresi,  “Even-Order  Aberration  Cancellation  in  Quantum 
Interferometry”,  Physical  Review  Letters,  v.  101,  233603  (2008). 

4)  While  investigating  the  physics  of  dispersion  cancellation  effects  we  have  formulated 
specific  conditions  enabling  total  spatial  dispersion  cancellation  to  all  orders.  The 
following  manuscript  has  been  submitted  for  publication: 


D.  S.  Simon  and  A.  V.  Sergienko  “Spatial-Dispersion  Cancellation  in  Quantum 
Interferometry”,  submitted  to  Physical  Review  A  (2009) 


As  a  final  step  in  our  research  efforts  on  this  project  we’ve  extended  the  notion  of  spatial 
dispersion  (aberration)  cancellation  to  the  case  of  “ghost  imaging”.  This  novel  imaging 
configuration  has  been  proposed  and  theoretically  analyzed  in  a  recent  manuscript  (D.  S. 
Simon  and  A.  V.  Sergienko  “Ghost  Imaging  with  Aberration-Cancellation”)  and  could 
open  the  way  for  enhanced-quality  aberration-free  imaging.  The  paper  been  submitted  for 
publication  in  Physical  Review  Letters. 

Conclusion: 

We  have  successfully  completed  the  original  research  program.  A  new  research  effect 
“Ghost  Imaging  with  Aberration-Cancellation”  has  been  discovered  in  the  process, 
leading  to  the  promise  of  significant  enhancements  in  the  fields  of  optical  imaging  and 
microscopy,  as  well  as  potential  improvements  to  phase-sensitive  optical  sensors. 

Our  research  has  led  to  five  research  papers  published  in  leading  refereed  physics 
journals,  with  one  additional  manuscript  recently  submitted  for  publication.  Nine 
conference  presentations  have  been  given,  including  two  Invited  Talks.  Four  Invited 
Lectures  have  also  been  given  at  leading  US  and  European  research  institutions,  and 
several  graduate  and  undergraduate  students  have  been  trained. 
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We  investigate  how  quantum  coincidence  interferometry  is  affected  by  a  controllable  manipulation  of  trans¬ 
verse  wave  vectors  in  type-II  parametric  down-conversion  using  adaptive  optics  techniques.  In  particular,  we 
discuss  the  possibility  of  spatial  walk-off  compensation  in  quantum  interferometry  and  an  effect  of  even-order 
spatial  aberration  cancellation. 
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I.  INTRODUCTION 

Quantum  entanglement  [1]  is  a  valuable  resource  in  many 
areas  of  quantum  optics  and  quantum  information  process¬ 
ing.  One  of  the  most  widespread  techniques  for  generating 
entangled  optical  states  is  spontaneous  parametric  down- 
conversion  (SPDC)  [2-5].  SPDC  is  a  second-order  nonlinear 
optical  process  in  which  a  pump  photon  is  split  into  a  pair  of 
new  photons  with  conservation  of  energy  and  momentum. 
The  phase-matching  relation  establishes  conditions  to  have 
efficient  energy  conversion  between  the  pump  and  the  down- 
converted  waves,  called  signal  and  idler.  This  condition  also 
sets  a  specific  relation  between  the  frequency  and  the  emis¬ 
sion  angle  of  down-converted  radiation.  In  other  words,  the 
quantum  state  emitted  in  the  SPDC  process  cannot  be  factor¬ 
ized  into  separate  frequency  and  wave-vector  components. 
This  leads  to  several  interesting  effects  where  the  manipula¬ 
tion  of  a  spatial  variable  affects  the  shape  of  the  polarization- 
temporal  interference  pattern.  For  example,  the  dependence 
of  polarization-temporal  interference  on  the  selection  of  col¬ 
lected  wave  vectors  was  studied  in  detail  in  [6]. 

Here  we  engineer  the  quantum  state  in  the  space  of  trans¬ 
verse  momentum  and  we  study  how  this  spatial  modulation 
is  transferred  to  the  polarization- spectral  domain  by  means 
of  quantum  interferometry.  We  will  focus  on  type-II  SPDC 
using  birefringent  phase  matching  since  the  correlations  be¬ 
tween  wave  vectors  and  spectrum  are  stronger  than  those 
employing  other  phase-matching  conditions. 

Our  aim  is  twofold.  From  one  point  of  view,  we  study  the 
effect  of  spatial  modulations  on  temporal  quantum  interfer¬ 
ence.  This  could  be  useful,  for  example,  in  quantum  optical 
coherence  tomography  (QOCT)  [7,8].  When  focusing  light 
on  a  sample  with  nonplanar  surface,  the  photons  will  acquire 
a  spatial  phase  distribution  in  the  far  field,  which  may  per¬ 
turb  the  shape  of  the  interference  dip.  Our  results  will  pro¬ 
vide  a  tool  to  understand  this  effect. 

From  a  second  point  of  view,  we  would  like  to  study  and 
characterize  spatial  modulation  as  a  tool  for  quantum  state 
engineering.  This  may  find  application  in  the  field  of  quan¬ 
tum  information  processing,  where  it  is  important  to  gain  a 
high  degree  of  control  over  the  production  of  quantum  en¬ 
tangled  states  entangled  in  one  or  more  degrees  of  freedom 
(hyperentanglement) . 

We  start  by  introducing  a  theoretical  model  of  a  type-II 
quantum  interferometer,  comprising  the  polarization,  spec¬ 


tral,  and  spatial  degrees  of  freedom  (Sec.  II).  A  modulation  in 
the  wave-vector  space  is  provided  by  an  adaptive  optical 
setup  and  equations  for  the  polarization-temporal  interfer¬ 
ence  pattern  in  the  coincidence  rate  are  derived.  In  Sec.  Ill, 
we  introduce  a  numerical  approach  for  practical  evaluation 
of  the  results  of  the  theoretical  model,  discussing  a  few  ex¬ 
amples  for  general  spatial  aberrations. 

In  Sec.  IV  we  will  highlight  and  discuss  theoretically  two 
interesting  special  cases.  The  first  one  is  the  possibility  of 
restoring  high  visibility  in  type-II  quantum  interference  with 
large  collection  apertures.  In  some  situations,  to  collect  a 
higher  photon  flux  or  a  broader  photon  bandwidth,  it  can  be 
useful  to  enlarge  the  collection  apertures  of  the  optical  sys¬ 
tem.  But  when  dealing  with  type-II  SPDC  in  birefringent 
crystals,  for  large  collection  apertures  the  effect  of  spatial 
walkoff  introduces  distinguishability  between  the  photons, 
leading  to  a  reduced  visibility  of  temporal  and  polarization 
quantum  interference.  We  will  show  that  high  visibility  can 
be  restored  with  a  linear  phase  shift  along  the  vertical  axis. 

The  second  effect  is  the  spatial  counterpart  of  spectral 
dispersion  cancellation  [9,10].  In  the  limit  of  large  detection 
apertures,  the  correlations  between  the  photons’  momenta 
will  cancel  out  the  effects  of  even-order  aberrations,  exactly 
as  in  the  limit  of  slow  detectors  the  frequency  correlations 
cancel  out  the  even-order  terms  of  spectral  dispersion.  The 
experimental  demonstration  of  this  effect  has  been  reported 
recently  [11]. 

As  we  proceed  from  the  general  case  of  Secs.  II  and  III 
into  the  specific  examples  of  Sec.  IV,  we  will  gradually  see 
that  optical  aberration  is  a  subject  with  two  very  different 
faces.  On  one  hand,  aberrations  in  optical  components  are 
normally  seen  as  undesirable,  since  they  lead  to  distortions  in 
imaging.  We  will  see  that  these  unwanted  spatial  modula¬ 
tions  may  to  some  extent  be  canceled.  On  the  other  hand,  we 
will  find  that  such  spatial  modulations  may  also  be  turned 
into  a  useful  tool:  by  deliberately  introducing  spatial  modu¬ 
lations  (in  effect,  purposely  adding  aberrations  in  a  con¬ 
trolled  manner),  we  can  produce  useful  effects  such  as  the 
restoration  of  visibility  mentioned  above.  The  examples  we 
provide  in  Sec.  IV  will  illustrate  these  two  aspects  and  will 
show  that  both  can  benefit  from  more  detailed  study  of  the 
interplay  between  spatial  modulations  and  spectral  correla¬ 
tions. 
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EIG.  1.  (Color  online)  Scheme  of  the  proposed  setup.  Horizon¬ 
tally  polarized  photons  from  type-II  SPDC  are  assigned  a  phase 
dependent  on  the  photon  transverse  momentum  (poiqo),  while  the 
vertically  polarized  ones  are  assigned  a  phase  4>e{qe)-  The  modu¬ 
lated  photons  enter  a  type-II  quantum  interferometer,  which  records 
the  coincidence  count  rate  as  a  function  of  the  delay  r  between  the 
photons  given  by  an  appropriate  delay  line. 

II.  THEORETICAL  MODEL 

Consider  the  scheme  in  Fig.  1.  A  laser  beam  pumps  a 
nonlinear  material  phase  matched  for  type-II  parametric 
down-conversion,  creating  a  pair  of  entangled  photons.  Each 
of  the  generated  photons  passes  through  a  Fourier-transform 
system,  and  then  enters  a  modulation  system  which  trans¬ 
forms  transverse  wave  vectors  for  the  horizontally  (H)  polar¬ 
ized  photon  according  to  the  transfer  function  Gi(qi)  and  for 
the  vertically  (V)  polarized  photon  according  to  G2(q2).  Af¬ 
ter  being  modulated  in  q  space,  the  photons  enter  a  type-II 
interferometer.  A  nonpolarizing  beam  splitter  (BS)  creates 
polarization  entanglement  from  the  polarization-correlated 
pair  emitted  by  the  source.  The  beams  at  the  output  ports  of 
the  beam  splitter  are  directed  toward  two  single-photon  de¬ 
tectors.  Two  polarizers  at  45°  before  the  BS  restore  indistin- 
guishability  in  the  polarization  degree  of  freedom.  An  adjust¬ 
able  delay  line  r  is  scanned  and  the  coincidence  rate  R{t) 
between  the  detection  events  of  the  two  detectors  is  recorded. 
An  aperture  is  placed  before  the  beam  splitter  to  select  an 
appropriate  collection  angle. 

A.  Notation 

Consider  a  monochromatic  plane  wave  of  complex  ampli¬ 
tude  £’(r)=£’o^“^*^*’,  with  r=(v,y,z).  For  a  given  wavelength 
X,  corresponding  to  a  frequency  H,  the  wave  vector  can  be 
split  into  a  transverse  component  ^  =  {qx^qy)  and  a  longitudi¬ 
nal  component  )S(q;H): 


k  =  [q.yS(q;a)]. 

(1) 

The  wave  number  is 

c 

(2) 

The  longitudinal  component  of  the  wave  vector  is 

yS(q,n)  =  -  |qp. 

(3) 

Therefore  the  electric  field  at  the  position  r  and  time  t 
written  as 

can  be 

E{r;t)  =  j  dqj  (/ll£(q,ri)e-'‘'''e-'^)^e'^, 

(4) 

where  p={x,y). 

In  paraxial  approximation  so  that 

,5) 

For  a  quasimonochromatic  wave  packet  centered  around  the 
frequency  Hq,  one  can  write  with  This 

expression  can  be  approximated  by 

I  P 

yS(q,a)-fco+--!^,  (6) 

Mo  2^0 

where  ko=k(ilo)  and  mo=(  is  the  group  velocity 

for  the  propagation  of  the  wave  packet  through  the  material. 


B.  State  generation 

Using  first-order  time-dependent  perturbation  theory,  the 
two-photon  state  at  the  output  of  the  nonlinear  crystal  can  be 
calculated  as 


dtH,{t)\0},  (7) 

where  the  interaction  Hamiltonian  is 

H,it)  =  f  J  drx^^\r)E^;\rM;\r,t)E^\r,t) •  (8) 

The  strong,  undepleted  pump  beam  can  be  treated  classically. 
Assuming  a  monochromatic  plane  wave  propagating  along 
the  z  direction, 

=  (9) 

The  signal  and  idler  photons  are  described  by  the  following 

quantum  field  operators: 

E]r\r,t)=  I  d(\j  \  wy) ,  (10) 


where  j=e,o. 

The  biphoton  quantum  state  at  the  output  plane  of  the 
nonlinear  crystal  is  [12] 


dv^iq,  v)dl(q,^o  +  v)dl{-  q,^o  -  *")|0)- 


(11) 


Two  photons  are  emitted  from  the  nonlinear  crystal,  one 
horizontally  polarized  (ordinary  photon)  and  the  other  verti¬ 
cally  polarized  (extraordinary  photon),  with  anticorrelated 
frequencies  and  emission  directions. 

In  the  case  of  a  single  bulk  crystal  of  thickness  L  and 
constant  nonlinearity  Xo^  the  probability  amplitude  for  having 
the  signal  photon  in  the  mode  (q,Ho+  and  the  idler  in  the 
mode  (-q,Ho-r')  is 


0(q,  p)  =  sine 


LA(q,  p) 
2 


^i[A{q,v)L/2] 


(12) 


For  type-II  collinear  degenerate  phase  matching,  the  phase- 
mismatch  function  A(q,  r')  can  be  approximated  to  be 
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A(q,v)  =  - v£)-l-Me2-q+-— ,  (13) 

kp 

where  D  is  the  difference  between  the  inverse  of  the  group 
velocities  of  the  ordinary  and  extraordinary  photons  inside 
the  birefringent  crystal  and  the  quadratic  term  in  q  is  due  to 
diffraction  in  paraxial  approximation.  The  remaining  term  is 
the  first-order  approximation  for  the  spatial  walkoff. 


C.  Propagation 

Consider  a  photon  described  by  the  operator  dj{q ,  H)  (po¬ 
larization  j=e,o,  frequency  H,  and  transverse  momentum  q). 
Its  propagation  through  an  optical  system  to  a  point  on  the 
output  plane  is  described  by  the  optical  transfer  function 
ify(x^,q;H).  In  our  setup,  the  field  at  the  detector  will  be  a 
superposition  of  contributions  from  the  ordinary  and  extraor¬ 
dinary  photons.  The  quantized  electric  fields  at  the  detector 
planes  are 

=  J  dqj  i(we'""^[//<,(x^,q;w)4(q,w) 

+  //o(x^,q;w)a„(q,w)], 


J  JqJ  d(oe^^^B[H^(xg,q;(o)dg(q,a)) 

+  q ;  o))d^iq,  w)] .  (14) 

The  probability  amplitude  for  detecting  a  photon  pair  at  the 
detector  planes,  with  space-time  coordinates  (x^Ja) 
(xg,^^),  is 


1.  State  engineering  section 

In  the  state  engineering  section,  each  of  the  two  branches 
consists  of  a  pair  of  achromatic  Fourier-transform  systems 
coupled  by  a  spatial  light  modulator  or  a  deformable  mirror. 
Each  Fourier-transform  system  consists  of  a  single  lens  of 
focal  length  /,  separated  from  the  optical  elements  before 
and  after  it  by  a  distance  /.  The  first  Fourier  system  maps 
each  incident  transverse  wave  vector  q  on  the  plane  to  a 
point  x(q)  on  the  Fourier  plane  11^: 

x(q)  =  -f-q,  kQ=—,  (17) 

ko  C 

where  /  is  the  focal  length  of  the  Fourier-transform  system. 
Since  we  assume  that  the  system  is  achromatic  for  a  certain 
bandwidth  around  a  central  frequency  Hq,  the  position  x(q) 
depends  only  on  q  and  not  on  o). 

The  spatial  modulator  assigns  a  different  amplitude  and 
phase  to  the  light  incident  on  each  point,  as  described  by  the 
function  G{x)  =  t{x)e^^^^\  Each  point  is  then  mapped  back  to 
a  wave  vector  on  the  plane  11^^^^  by  the  second  achromatic 
Fourier-transform  system. 

Using  the  formalism  of  Fourier  optics  [13],  the  transfer 
function  between  the  planes  11  ^out  be  calculated 
to  be 


The  corresponding  momentum-transfer  function  is 

7 


^^i(qi>q3)  =  G 


-qi 


^qi-qa)- 


(19) 


A(xA,XB-,tA^tB)  =  (0|£7(xa>1a)£b  ^(xb.^b)I*A)-  (15) 

For  the  biphoton  probability  amplitude  we  get 

A{XA,XB\tA,tB)^  /  dq^dq^dcj^dcj/^^q^,  q^ ; 

X  iH,{xA,  q^ ;  0)^)H„{xb,  q^ ;  w 
-I-  //o(x^,qo;  w„)//<,(xs,q,; 

(16) 

This  probability  amplitude  represents  the  superposition  of 
two  possible  events  leading  to  a  coincidence  count  in  the 
detectors: 

(1)  the  V  polarized  photon  with  momentum  q^  and  fre¬ 
quency  cOg  going  through  the  lower  branch  to  arrive  at  posi¬ 
tion  in  detector  A,  while  the  H  polarized  photon  with 
momentum  q^  and  frequency  goes  through  the  upper 
branch  to  arrive  at  position  in  detector  B\  and 

(2)  the  V  polarized  photon  with  momentum  q^  and  fre¬ 
quency  cOg  going  through  the  lower  branch  to  arrive  at  posi¬ 
tion  Xb  in  detector  B,  while  the  H  polarized  photon  with 
momentum  q^  and  frequency  goes  through  the  upper 
branch  to  arrive  at  position  in  detector  A. 

Since  the  superposition  is  coherent,  there  are  quantum- 
interference  effects  between  the  two  probability  amplitudes. 


2.  Interferometer 

After  the  plane  11^^^^  the  two  photons  enter  a  type-II  quan¬ 
tum  interferometer.  Each  propagates  in  free  space  to  a  bire¬ 
fringent  delay  line  and  a  detection  aperture  p{x)  to  be  finally 
focused  to  the  detection  planes  by  means  of  lenses  of  focal 
length  /q.  Following  the  derivation  in  [6]  the  transfer  func¬ 
tion  is 


H^ix, 


i,q;w)  =  j 


h(xi,Xi;co)e^^'^^dxi 


-  ^i(oj/c)(di+d2+fo) 


exp 


2c/o  \/o 


\c/o 


(20) 


where  P(q)  is  the  Fourier  transform  of  |/7(x)p. 

A  combination  of  the  two  different  stages  is  described  by 
the  transfer  function 


f 

—  C 

ko 


H2(x^,qa,(^), 


(21) 


where  the  two  functions  Gi(q)  and  G2(q)  are  the 
momentum-transfer  functions  which  describe  the  modulation 
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imparted,  respectively,  on  the  ordinary  and  the  extraordinary 
photons. 


D.  Detection 

Since  the  single-photon  detectors  used  in  quantum  optics 
experiments  are  slow  with  respect  to  the  temporal  coherence 
of  the  photons  and  their  area  is  larger  than  the  spot  onto 
which  the  photons  are  focused  by  the  collection  lens,  we 
integrate  over  the  spatial  and  temporal  coordinates.  There¬ 
fore  the  coincidence  count  rate  expressed  in  terms  of  the 
biphoton  probability  amplitude  is 


R{t)^  j  dx^  j  dxjij  dtA  j  dtjilA(xA,Xji:tAdB)l^.  (22) 


Following  the  derivation  described  in  Appendix  A,  one 
gets 


R(t)=Ro 


1-A 


Wm(t) 


where  A(x)  is  the  triangular  function 


A(v)  = 


1  -  |v| 

0, 


kl  ^  1 

|v|  >  1. 


(23) 


(24) 


Therefore,  the  coincidence  count  rate  R{t)  is  given  by  the 
summation  of  a  background  level  Rq  and  an  interference  pat¬ 
tern  given  by  the  triangular  dip  A(1  -^)  that  one  gets  when 
working  with  narrow  apertures,  modulated  by  the  function 
1T^(t)  which  depends  on  the  details  of  the  adaptive  optical 
system. 

The  expressions  for  Rq  and  W^(t)  are 


EIG.  2.  (Color  online)  Example  of  the  numerical  approach 
adopted  to  evaluate  Eqs.  (25)  and  (26).  The  spatial  modulation  sur¬ 
face  is  discretized  in  sufficiently  small  squares  over  which  the  phase 
is  averaged. 

III.  NUMERICAL  SOLUTIONS  FOR  A  GENERAL 
PHASE  SHIFT 


Rn 


sinc[MLe2  •  (q  -  q')] 


Zq' j^-4^L/2)e24q-q0^42^/i/y[|qp-|q'P] 
XP^iq  -  q')PB{-  q  +  q')  (25) 


and 


MLe2-(q  +  q')A 


I 

XI-  — 
DL 


Gllfq  G,(fq'  G;  -fq 


ko 


ko 


X  G- 


_  L-/(M/D)7«2-(q-q')^42^/i/y[|qp-|q'p] 

ko^  r 


x-P/ik  +  q'M-Cq  +  q')]-  (26) 

In  the  following  we  will  assume  there  is  spatial  modula¬ 
tion  only  on  one  of  the  photons;  therefore  we  set  G2(q)  =  1. 


Numerically  solving  for  the  quantities  in  Eqs.  (25)  and 
(26)  in  the  case  of  a  general  aberration  may  be  computation¬ 
ally  demanding.  Here,  we  propose  an  approximation,  valid  in 
the  case  where  the  function  G(x)  changes  smoothly  over  the 
mirror  surface,  as  it  is  the  case  in  experimentally  relevant 
situations.  This  model  is  also  interesting  from  the  practical 
point  of  view,  since  in  many  cases  adaptive  optical  systems 
are  implemented  using  spatial  light  modulators  or  segmented 
deformable  mirrors,  where  the  modulation  surface  is  parti¬ 
tioned  into  small  pixels. 

Suppose  we  partition  the  Fourier  plane  H^  into  small 
squares  (pixels)  of  side  d  (Fig.  2).  Let  us  define  the  rectan¬ 
gular  function 


n(x)  = 


if  |x|  >  I 

if  |x|  < 


(27) 


The  pixel  {l,m)  is  identified  by 


=  n 


X 

n 

y 

~  1 

d 

id  J 

(28) 


selecting  the  area 


{l-j)d<x<{l+^)d, 
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(m  -  <  y  <  (m  +  |)j.  (29) 

We  approximate  the  value  of  the  phase  in  each  square  by 
the  mean  value  of  the  actual  phase  within  the  square: 


y 

“  +  / 

n 

Id  \ 

[d  J 

That  is, 


a  dxdy(p(x,y)n 


^i(p(x,y)  _  ^  ^i(pi  ,^n[(x/d)+l]n[(y/d)+m]^ 
l,m 


(31) 


In  this  case  (see  Appendix  B  for  a  justification) 
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2  ^i(f>l,m^[(x/d)+l]n[(y/d)+m]  _  ^ 
l,m  l,m 


y 

“  +  / 

n 

-  +  m 

[d  J 

d 

(32) 


(30)  Substituting  this  expression  into  Eq.  (23)  and  collecting  the 
integrations,  one  finds 


R(t)  ~  2  2  (33) 

l,m  k,/uL 


where 


=  j  dq^j  dQJl 


f  j 


n 


/ 

—Qx-^ 


gmdilkp)(ql-Ql)P[q^+QM-(‘lx+Qx)'\ 


and 


4,m('^) 


k/ 


XSinc 

Performing  the  integrations  one  gets 

=  J  dxP(x)A 

and 


=  j  dqyj  dQyll 

ML(qy  +  Qy)A{  1 


n 

2t 


f  n 


DL 


P[.Qy  +  Gy]^[“  (^y  +  Gy)] 


fx  .  , 


j  2ddi 

sinci - xA 

I  / 


fx  .  . 


J{ddilf){l-\)x 


/ 


Im/x('^)=  I  dxP(x)A 


fx 

kd 


-  (m  +  /x) 


sine] 

1  MLxA 

2t 

1-  — 

\ 

J  1 

[ 

DL. 

J 

\lkdlldy  M  . 

X  sinci -  — V-  — T  A 

i  /  \  D  ' 


—  -  (m  +  /x) 
kd 


^ i {kdlf) \_{2d^ Ik  )x- {MID)  r] 


A  similar  expression  can  be  found  for  the  background  coincidence  rate: 


l,m  \,fji 


where 


and 


=  j  dxP(x)A 


fx  .  . 


(  2dd^ 

sinci  - xA 

I  / 


fx  .  . 


Rml^j  dxP(x)A 


—  -{m-  fj) 
kd  ^ 


\  2ddi 

smc(MLv)smc|  ^  xA 


—  -{m-  ff) 
kd 


i{ddilf){l+\)x 


i[{ddilf){m+ fji)-MLI2]x 


(34) 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


The  advantage  of  our  numerical  approach  is  that  one  can 
calculate  and  tabulate  the  functions  and  /^^(r) 

for  a  given  configuration,  determined  by  the  focal  length  /, 
the  shape  of  the  detection  apertures,  the  width  of  the  deform¬ 
able  optics,  and  the  distance  between  the  crystal  and  the 


detectors.  Then,  to  calculate  the  shape  of  the  interference 
pattern  for  a  certain  phase  distribution  on  the  adaptive  optics, 
one  just  needs  to  change  the  coefficient  of  a  linear  combina¬ 
tion  of  the  tabulated  functions.  This  can  be  a  helpful  tool  for 
studying  the  effect  of  specific  aberrations  on  the  temporal 
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EIG.  3.  (Color  online)  Examples  of  the  shape  of  the  polarization 
quantum-interference  dip  with  two  different  spatial  phase  modula¬ 
tions  (in  black  the  unperturbed  dip,  in  red  the  modulated  one).  In 
the  upper  figure  the  effect  of  a  small  amount  of  coma  along  the 
vertical  axis  is  shown.  In  the  lower  figure  a  more  complicated  su¬ 
perposition  of  aberrations  affects  dramatically  the  shape  of  the  dip. 

interference  or  to  engineer  the  shape  of  the  Hong-Ou-Mandel 
(HOM)  dip. 

Some  examples  of  how  the  temporal  quantum- 
interference  dip  is  modulated  by  a  generic  spatial  phase  shift 
are  reported  in  Fig.  3,  for  coma  (upper  plot)  and  a  superpo¬ 
sition  of  several  different  aberrations  (lower  plot).  The  inter¬ 
ference  visibility  clearly  degrades  in  presence  of  wave-front 
aberrations. 


IV.  PARTICULAR  CASES 

In  this  section  we  will  discuss  the  quantum-interference 
pattern,  described  by  Eq.  (23),  for  a  few  simple  cases.  First 
we  will  consider  the  case  when  no  spatial  modulation  is  as¬ 
signed  to  the  photons  and  Eq.  (23)  will  reduce  to  the  results 
already  described  in  the  literature  for  quantum  interferometry 
with  multiparametric  entangled  states  from  type-II  down- 
conversion  [6].  Then  we  will  examine  the  effect  of  a  linear 
phase,  describing  its  implications  for  the  compensation  of  the 
spatial  walkoff  between  the  two  photons.  Finally  we  will 
describe  what  happens  in  the  approximation  of  sufficiently 
large  detection  apertures,  introducing  the  effect  of  even-order 
aberration  cancellation. 

A.  No  phase  modulation 

Applying  no  phase  modulation,  our  equations  reduce  to 
the  ones  derived  in  [6].  Particularly  we  find 


7?o  =  -Pa(0)Pb(0) 


(41) 


and 


Wg{t)  =  Sine 


M^Lk„  /  2t 

- ^tA  1  -  — 

2diD  \  DL 


Mk^ 

Mk„ 

X 

- ^^2 

2diD 

Pb 

- - ^^2 

2diD 

(42) 


The  shape  of  the  interference  pattern  is  essentially  given  by 
the  product  of  the  triangular  function  by  the  Fourier  trans¬ 
form  of  the  aperture  function,  centered  at  t=0.  For  physi¬ 
cally  relevant  parameters  the  sine  function  is  almost  flat  in 
the  region  where  the  triangular  function  is  not  zero. 

To  get  an  analytic  result  one  may  assume  Gaussian  detec¬ 
tion  apertures  of  radius  Rq  centered  along  the  system’s  opti¬ 
cal  axis: 


p(x)  =  e 

In  this  case  the  solution  is  quite  simple: 


R{t)=R, 


(43) 


(44) 


with 


Ti  = - (45) 

kpMRc 

Typically,  sharp  circular  apertures  are  used  in  experi¬ 
ments.  In  this  case,  the  function  P[q]  is  described  in  terms  of 
the  Bessel  function  /^(v).  For  a  circular  aperture  of  radius  R, 


/^[q]  = 


■/i(2^|q|) 

^|q| 


(46) 


However  the  Gaussian  approximation  is  still  a  good  one  if 
the  width  Rq  of  the  Gaussian  is  taken  to  roughly  fit  the 
Bessel  function  (of  width  R)\  in  our  case  we  take  Rq 
=R/(2^). 

Therefore  Eq.  (44)  is  still  approximately  valid  in  the  case 
of  sharp  circular  apertures,  just  taking 


4V2t/i£> 

kpMRg 


(47) 


Mathematically,  in  Eq.  (44)  the  interference  pattern  is 
given  by  the  multiplication  of  a  triangular  function  centered 
at  t=DLI2  by  a  Gaussian  function  centered  at  t=0.  The 
width  of  the  Gaussian  function  decreases  with  increasing 
radius  of  the  aperture  Rp.  Therefore,  in  the  small- aperture 
approximation,  the  width  of  the  Gaussian  is  so  large  that  it  is 
approximately  constant  between  t=0  and  t=DLI2,  giving 
the  typical  triangular  dip  found  in  quantum-interference  ex¬ 
periments.  On  the  other  hand,  increasing  the  aperture  size, 
the  width  of  the  Gaussian  function  decreases,  reducing  the 
dip  visibility  (see  Fig.  4).  Physically,  this  can  be  explained 
by  the  fact  that  by  increasing  the  aperture  size  we  let  more 
wave  vectors  into  the  system,  and  so  the  spatial  walkoff  in 
type-II  interferometry  introduces  distinguishability,  reducing 
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EIG.  4.  On  the  right  side  we  can  see  the  interference  patterns 
with  three  different  detector  aperture  sizes:  the  corresponding  aper¬ 
ture  functions  are  shown  on  the  left  side. 

the  interference  visibility.  Enlarging  the  aperture  sizes  is  of¬ 
ten  useful  in  practice,  for  example,  to  get  a  higher  photon 
flux.  Moreover,  since  in  the  SPDC  process  different  fre¬ 
quency  bands  are  emitted  at  different  angles,  it  may  be  nec¬ 
essary  to  open  the  detection  aperture  in  applications  where  a 
broader  bandwidth  is  needed.  This  is  clearly  a  problem  when 
using  type-II  phase  matching  in  birefringent  crystals,  since 
the  visibility  of  temporal  and  polarization  interference  gets 
drastically  reduced. 


B.  Linear  phase  shift 


EIG.  5.  Quantum-interference  pattern  for  different  detector  ap¬ 
erture  sizes,  introducing  a  linear  modulation  of  the  deformable  mir¬ 
ror,  in  order  to  restore  the  indistinguishability  between  the  photons, 
decreased  by  the  spatial  walkoff  in  the  generation  process. 


koML 
2/  ' 


(52) 


In  the  case  of  a  reflective  system,  in  which  the  phase 
modulation  is  implemented  by  means  of  a  deformable  mirror 
(Fig.  5),  tilted  by  an  angle  6, 


<p(x)  =  2kQ  tan  0y-Siyy.  (53) 

Therefore,  the  amount  of  tilt  necessary  to  restore  high  vis¬ 
ibility  is 


Suppose  now  we  introduce  a  linear  phase  function  with 
the  spatial  light  modulator,  along  the  direction  Sj, 


^(x)  =  Si-x, 


(48) 


we  get 


M^Lk„  1 

2t\ 

Sine 

- ^tA\ 

1-  — 

2d^D  \ 

V  DL 

Mk„ 

Mk„  ^ 

X 

Pb 

(49) 


If  we  compare  Eq.  (49)  with  Eq.  (44)  we  can  see  that  the 
structure  is  the  same.  We  again  have  a  triangular  function 
centered  at  t=DLI2,  along  with  two  aperture  functions.  But 
this  time,  instead  of  being  centered  at  t=0,  the  aperture  func¬ 
tions  can  be  shifted  at  will  along  the  r  axis.  Suppose  we  now 
apply  a  tilt  along  the  y  axis  {si^=0).  The  modulation  function 
is  then  shifted  to 


'^center  j  (50) 

KqM 

To  get  the  highest  possible  visibility,  the  center  of  the  modu¬ 
lation  function  must  be  matched  to  the  center  of  the  triangu¬ 
lar  dip, 


^center 


DL 

T’ 


(51) 


SO  that 


tan  - .  (54) 

4/ 

In  the  case  of  a  1.5  mm  crystal,  with  M= 0.0723  (pump  at 
405  nm,  SPDC  at  810  nm)  and  lenses  with  focal  length  of  20 
cm  in  the  4/  system,  we  get 

^=0.14  mrad.  (55) 


C.  Large-aperture  approximation 

If  the  detection  apertures  are  large  enough  for  the  func¬ 
tion  to  be  successfully  approximated  by  a  delta  function,  we 
get 


— n  L-«(2M/D)7€2  q 

Cr 


(56) 


Suppose  that  the  spatial  modulator  is  a  circular  aperture 
with  radius  r,  with  unit  transmission  and  phase  modulation 
described  by  the  function  ^(x). 


0  if  |x|  >  r 
1^1  ^  ^ 


(57) 


In  this  case  the  function  (p(x)  can  be  expanded  on  a  set  of 
polynomials  which  are  orthogonal  on  the  unit  circle,  such  as 
the  Zernicke  polynomials: 


=  22  <P„mRnip)^Os{me), 

n  m 
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FIG.  6.  Plot  of  the  values  of  \a,J^,  for  /,\=-100, ...  ,+100. 
The  radius  of  the  detection  apertures  is  R  =  5  mm,  the  distance  be¬ 
tween  the  exit  plane  of  the  modulation  section  and  the  detection 
apertures  is  d- 1  =  l  m,  and  the  size  of  the  modulation  pixels  is  d 
=  15  yum.  Clearly  only  the  diagonal  elements  are  nonzero,  i.e.,  the 
ones  for  which  \  =  -l.  In  this  situation  the  effect  of  even-order  ab¬ 
erration  cancellation  is  present. 


m  =  -n,-n  +  2,-n-\-A,  ...  (58) 

where  q=(p  cos  ^,p  sin  6).  To  calculate  ^(-q)  we  note  that 
-q  =  [p  cos(^+7r),p  sin(^-i-7r)],  so 

<p(-q)  =  ESC(p)cosM0+77)].  (59) 


If  m  is  even  then  cos[m(^+7r)]  =  cos(m^);  otherwise  if  m  is 
odd  cos[m(^+7r)]=  -cos  mO).  Therefore 

^)(q)  -  «p(- q)  =  2^  2  ?>«m^;!'(p)cos(mP).  (60) 

n  m  odd 

So,  only  the  Zernicke  polynomials  with  m  odd  contribute 
to  the  shape  of  the  interference  pattern.  This  effect  is  the 
spatial  counterpart  of  the  dispersion  cancellation  effect,  in 
which  only  the  odd-order  terms  in  the  Taylor  expansion  of 
the  spectral  phase  survive.  The  experimental  demonstration 
of  this  effect  was  recently  reported  in  [11]. 

An  interesting  question  is  how  large  the  detection  aper¬ 
tures  should  effectively  be,  in  order  to  obtain  the  even-order 
aberration  cancellation  effect.  According  to  the  numerical  ap¬ 
proach  proposed  in  Sec.  IV,  the  even-order  aberration  can¬ 
cellation  effect  manifests  itself  in  the  limit  where  P(v) 
~  ^(v),  so  that 

+  (61) 

In  Fig.  6,  a  plot  of  the  value  for  a/  ^  is  shown  for  typical 
values  of  the  relevant  experimental  parameters  (detection  ap¬ 
erture  radius  R  =  5  nnn,  detection  distance  di  =  \  m,  and  size 
J=0.1  mm  of  each  pixel  in  the  Fourier  plane  of  the  adaptive 
optical  system).  Clearly,  only  the  diagonal  elements  (the  ones 
for  which  l=-X)  are  significant,  suggesting  that  the  effect  of 
even-order  aberration  cancellation  may  be  observable  for 
most  typical  experimental  parameters. 

To  get  an  idea  of  what  happens  for  different  experimental 
conditions,  we  can  compute  the  ratio  between  the  intensities 


d  =  0.5  mm 


FIG.  7.  (Color  online)  Plot  of  the  ratio  po  between  the  intensities 
of  the  nondiagonal  coefficient  a^i  and  the  diagonal  coefficient 
as  a  function  of  the  radius  of  Gaussian  detection  apertures,  R,  for 
different  values  of  the  distance  between  the  exit  plane  of  the  modu¬ 
lation  section  and  the  detection  apertures  (Ji  =  1,10, 100  cm).  On 
the  upper  plot  the  size  of  the  modulation  pixel  d  is  J=0.5  mm, 
while  in  the  second  case  it  is  d=Q.15  mm.  Clearly,  for  experimen¬ 
tally  interesting  cases,  the  off-diagonal  coefficient  is  at  least  3  or¬ 
ders  of  magnitude  smaller  than  the  diagonal  one,  leading  to  even- 
order  aberration  cancellation. 


of  the  nondiagonal  coefficient  and  the  diagonal  coeffi¬ 
cient 


i».iP 


(62) 


The  lower  the  value  for  po  is,  the  less  significant  the  coeffi¬ 
cients  for  X  A  -/  are:  the  even-order  aberration  cancellation 
effect  will  therefore  manifest  itself  more  clearly. 

Values  for  po  are  shown  in  Fig.  7  for  two  different  cases. 
In  both  pictures,  the  value  of  po  is  shown  as  a  function  of  the 
Gaussian  detection  aperture  radius  R,  for  three  different  val¬ 
ues  of  the  distance  between  the  plane  113  the  detection 
lenses  Jj.  In  the  upper  panel,  the  size  of  each  small  square  in 
which  the  spatial  phase  is  assumed  to  be  constant  is  d 
=0.5  mm,  while  in  the  second  case  it  is  J=0.25  mm.  In 
both  cases  po  is  significantly  smaller  than  1,  and  it  becomes 
smaller  and  smaller,  increasing  the  value  of  the  detection 
aperture  radius.  However,  po  is  smaller  for  larger  values  of  J, 
implying  that  the  spatial  variability  of  the  modulation  phase 
plays  a  role  in  the  degree  of  even-order  cancellation  of  the 
modulation  itself. 

It  turns  out  that  for  the  aberration  cancellation  effect  to 
appear,  it  is  in  fact  only  necessary  for  one  aperture  to  be 
large  and  for  one  detector  to  be  integrated  over.  This  is  suf¬ 
ficient  to  produce  the  transverse-momentum  delta  functions 
that  lead  to  even-order  cancellation.  To  demonstrate  this,  we 
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can,  for  example,  consider  the  case  where  the  aperture  at  B  is 
large,  and  the  detector  at  B  is  integrated  over,  while  the  ap¬ 
erture  at  A  is  taken  to  be  finite,  with  detector  A  treated  as 
pointlike.  The  location  of  the  pointlike  detector  will  hence¬ 
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forth  simply  be  denoted  as  x,  and  we  continue  to  work  within 
the  quasimonochromatic  approximation.  If  we  integrate  only 
over  Vg,  leaving  v  unintegrated,  then  it  is  straightforward  to 
show  that  the  analogs  of  Eqs.  (A3)  and  (A4)  are 


(63) 


,v)^  q'  +  ^\pb(-  q'  -  q)  +  q(-  q+^\Q'(+  q'  +  ^\pB{q-  q')\ . 


cfo 


cfo 


cfo/  \  cfj 


(64) 


We  have  defined  Q  and  Q'  to  be  the  Fourier  transforms, 
respectively,  of  and  /7^.  is,  as  before,  the  Fourier  trans¬ 
form  of  \pb\^-  We  now  let  aperture  B  become  large,  so  that 

the  function  Pg  goes  over  to  a  delta  function.  For  Gi{x) 
=  G2(v)  =  1,  we  can  substitute  these  results  into  the 

coincidence  rate  (which  will  now  be  a  function  of  both  r  and 
the  position  v  of  detector  A),  and  carry  out  the  q'  and  v 
integrals.  For  the  modulation  term,  we  find 

7?^(x,  t)  =  P(x,  t)  -  Roix)  (65) 


+  2  “9^  + 


cfo  / 


cfo 


(66) 


Here,  we  have  used  the  fact  that  the  Fourier  transform  of 
/7^(v)  equals  the  complex  conjugate  of  the  Fourier  transform 

of  /7^(-v),  in  order  to  write  Q'  in  terms  of  Q.  We  see  from 
the  presence  of  the  factor  that  even-order  aber¬ 

ration  cancellation  occurs  even  though  one  aperture  is  finite 
and  the  corresponding  detector  is  pointlike.  This  point  may 
be  of  importance  in  future  attempts  to  produce  aberration- 
canceled  imaging. 


V.  CONCLUSIONS 

Summarizing,  in  this  paper  we  have  carried  out  a  theoret¬ 
ical  study  of  the  relation  between  the  wave-front  modulation 
of  the  entangled  SPDC  photons  and  the  shape  of  the  result¬ 
ing  temporal  quantum-interference  pattern.  Due  to  the  mul- 
tiparametric  nature  of  the  generated  entangled  states,  the 
modulation  on  the  spatial  degree  of  freedom  can  affect  the 
shape  of  the  polarization-temporal  interference  pattern  in  the 
coincidence  rate.  Our  aim  is  twofold:  from  one  side  we  want 


to  study  the  effect  of  wave-front  aberration  on  quantum  in¬ 
terferometry,  and  from  the  other  we  want  to  discuss  a  way  to 
engineer  multiparametrically  entangled  states. 

We  have  introduced  a  theoretical  model  for  calculation  of 
the  shape  of  the  polarization-temporal  interference  pattern 
given  a  certain  general  phase  modulation  in  the  crystal  far 
field,  assuming  as  a  free  parameter  the  shape  and  the  dimen¬ 
sion  of  the  collection  apertures.  Using  a  numerical  method  to 
study  the  resulting  equation  has  shown  that  for  typical  ex¬ 
perimental  cases  the  hypothesis  of  large  apertures  can  be 
assumed  to  be  valid.  In  such  an  approximation,  only  the  odd 
part  of  the  assigned  phase  modulation  affects  the  shape  of  the 
interference  pattern.  This  effect  has  recently  been  demon¬ 
strated  experimentally  [11]. 

Moreover,  it  is  often  useful  in  experiments  to  enlarge  the 
collection  aperture  in  order  to  collect  a  higher  photon  flux 
and  larger  optical  bandwidth.  But  when  working  with  type-II 
birefringently  phase-matched  down-conversion,  spatial 
walkoff  between  the  emitted  photons  introduces  distinguish- 
ability  between  the  two  possible  events  that  can  lead  to  co¬ 
incidence  detection,  reducing  the  visibility  of  quantum  inter¬ 
ference.  Such  walkoff  can  be  compensated  for  with  a  linear 
phase  shift  in  the  vertical  direction,  restoring  high  visibility. 
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APPENDIX  A:  SKETCH  OF  DERIVATION  OF  EQ.  (23) 


In  this  appendix  we  sketch  the  major  steps  for  the  derivation  of  Eq.  (23).  Substituting  Eq.  (21)  into  Eq.  (16)  and  the  result 
into  Eq.  (22),  one  finds  the  following  expressions  for  Rq  and  Wq{t): 

Ro  =  j  (Al) 

c?qfi;q'fi;i4)*(q,v)$(q',-v)G;(^q)Gi(^q')G;(-^q)G2(-^q')mq,q',i^),  (A2) 

where 


W<°)(q,q',v) 


I 


dx^(/xsH*(x^,q,v)H*(xs,-  q,-  v)H(x^,q',v)H(xs,- q',-v) 


+  H*(x^,-q,-  v)H*(xs,q,v)H(x^,-q',-  v)H(xs,q',v) 


(A3) 


and 


Wiq 


dxAdxBH*{xA,q,v)H*(xB-  q,-  v)H{xa,- q' ,v)H{xB,q'  -  v) 


f  r. 

'z  1 

dzU 

-h 

J 

.L  2 

+  H*{xj^-q,-  v)H*{xB,q,v)H{xA,q'  -  v)H{xB-q',v). 

The  angular  and  spectral  emission  function  <I)(q,  v)  is  given  by 

1  1 

^-iA{q,v)z 

Performing  the  integrals  over  the  spatial  coordinates  dx^  and  x^,  one  gets 

#'»(q,q',n)  =  V^[(q  _  q')]^^-  (q  -  q')]  +  Pa[-  (Q  -  q')]^B[(q  -  q')]} 

and 

lT(q,q',  n)  =  ]{p^[(q  +  q')]Pg[-  (q  +  q')]  +  (q  +  q')]PB[iq  +  q')]}- 


(A4) 

(A5) 

(A6) 

(A7) 


Einally,  use  of  the  integral  representation  for  the  sine  func¬ 
tion  [Eq.  (A5)]  allows  the  v  integration  to  be  carried  out,  but 
at  the  expense  of  introducing  two  integrations  over  a  pair  of 
new  parameters  (say,  z  and  z')-  Note  the  following  relation, 
which  can  easily  be  verified  by  sketching  the  functions  on 
the  left-hand  side: 


APPENDIX  B:  JUSTIFICATION  OF  EQ.  (32) 

Suppose  we  have  a  set  A,  which  can  be  partitioned  into  a 
collection  of  disjoint  subsets  A^,  with  k  =  1 ,2, . . .: 

UAy^  =  A,  A]^n  Ai  =  (f)  if  ki^l.  (Bl) 

k 


n[v]n[v  -  a]  = 


1 

1 

0 


if  -  l<a<0,  -|<x<|  +  a 

if0<a<l, 

otherwise. 


(A8) 


Erom  this,  it  follows  that 


To  each  set  we  can  associate  a  characteristic  function. 


Xkix)  = 


X  e 

X  f  Ai, 


(B2) 


such  that 


2  Xkix)  =  XAix),  Xkix)xiix)  =  Si^iXkix) ,  (B3) 

k 


J  n[v]n[v  -  a\dx  =  A(a) ,  (A9) 

where  A(a)  is  the  triangle  function.  These  facts  allow  us  to 
carry  out  the  two  z  integrations  that  arise  from  the  sine  func¬ 
tion,  leading  to  the  result  shown  in  Eq.  (23). 


where  Xa  is  the  characteristic  function  for  the  full  set. 


XAiA  = 


V  G  A 

V  $  A. 


(B4) 


The  term  assumes  the  value  of  for  Xk{^)  =  \  and 

the  value  of  1  for  ;^^(x)  =  0  [l-Xk(^)  =  ^^^  so 
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exp 


k 


n  =  n  {1[1  -  Xkix)^  +  e''^^Xkix)]  =  11  [1  +  DXkl 

k  k  k 


If  we  express  the  first  few  terms  we  get 


(B5) 


n  [1  +  -  Dxk]  =  [1  +  -  1)a:i][1  +  -  l);t2] ' ' ' 

k 

=  1  +  (e'^1  -  Dxi  +  (e'^2 -1)X2+---  +  +  (g'^i  _  i)(e'>3 _  D^^x^ 

+  . . .  +  (e''0i  _  i)(e'>i  _  i)(e«>i  _  l)xiX2X3  +  -  \)xiX2Xa  +  ' ' '  • 

So  that  in  the  end 


exp 


<t>kXkix) 

k 


=  i  +  E[(e‘^^-ikJ  =  i  +  2. 

k  k 


'^Xk-^Xk^^‘ 

k  k 


%• 


(B6) 

(B7) 


Since  the  square  sets  we  have  used  in  Sec.  IV  satisfy  Eq.  (Bl),  then  the  result  expressed  in  Eq.  (B5)  is  valid  for  our  case. 
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We  report  the  hrst  experimental  demonstration  of  even-order  aberration  cancellation  in  quantum 
interferometry.  The  effect  is  a  spatial  counterpart  of  the  spectral  group  velocity  dispersion  cancellation, 
which  is  associated  with  spectral  entanglement.  It  is  manifested  in  temporal  interferometry  by  virtue  of  the 
multiparameter  spatial-spectral  entanglement.  Spatially  entangled  photons,  generated  by  spontaneous 
parametric  down-conversion,  were  subjected  to  spatial  aberrations  introduced  by  a  deformable  mirror  that 
modulates  the  wave  front.  We  show  that  only  odd-order  spatial  aberrations  affect  the  quality  of  quantum 
interference. 


DOI:  10.1 103/PhysRevLett.  101. 233603 

The  nonlinear  optical  effect  of  spontaneous  parametric 
down-conversion  (SPDC)  has  been  a  reliable  source  of 
entangled-photon  states  for  the  last  30  years.  A  photon  of 
the  pump  radiation  has  a  random  chance  to  be  converted 
into  two  photons  with  lower  energy  while  traveling  though 
the  nonlinear  material.  Conservation  of  energy  and  mo¬ 
mentum  governs  the  spatial  and  spectral  state  of  the  down 
converted  light.  In  the  case  of  a  monochromatic  pump 
beam,  energy  conservation  leads  to  strong  anticorrelation 
between  the  frequency  components  of  signal  and  idler 
wave  packets.  This  symmetric  superposition  of  all  possible 
anticorrelated  frequencies  with  respect  to  the  degenerate 
frequency  of  signal  and  idler  waves  gives  rise  to  frequency 
entanglement. 

Even-order  dispersion  cancellation  is  among  the  most 
interesting  consequences  of  frequency  entanglement  [1,2]. 
If  one  of  the  two  photons  of  an  entangled  pair  travels 
through  a  dispersive  material  and  both  photons  are  com¬ 
bined  on  a  beam  splitter  in  a  Hong-Ou-Mandel  configura¬ 
tion  [3],  then  the  rate  of  coincidences  between  the  counts  of 
two  single-photon  detectors  placed  at  the  output  ports 
depends  on  the  odd-order  dispersion  terms  only  when 
observed  as  a  function  of  the  path  difference  between  the 
two  arms  before  the  beam  splitter.  The  detrimental  effect  of 
even-order  dispersion  (such  as  group  velocity  dispersion), 
which  leads  to  the  wave  packet  broadening,  is  canceled. 
This  has  been  exploited  in  several  applications  such  as  the 
measurement  of  photon  traveling  time  trough  a  material 
[4],  and  improving  the  accuracy  of  remote  clock  synchro¬ 
nization  [5].  Optical  coherence  tomography  [6,7]  has  also 
benefited  from  this  nonclassical  effect.  This  quantum  effect 
has  inspired  recent  developments  of  classical  nonlinear 
optical  systems  mimicking  dispersion  cancellation  [8,9]. 

The  wave  vector  of  a  monochromatic  wave  at  a  given 
frequency  ft  has  a  bidimensional  transverse  wave  vector  q 
(in  the  plane  orthogonal  to  the  propagation  direction)  and  a 

longitudinal  component  Ac(q,  ft)  =  V[^^(ft)ft^/<^^]  “  IqP- 
In  parametric  down-conversion  with  a  plane-wave  pump. 


PACS  numbers:  42.50.Dv,  03.67.Bg,  42.30.Kq,  42.50.St 

momentum  conservation  leads  to  anticorrelation  of  the 
transverse  wave  vector  components  [10].  This  analogy 
with  frequency  anticorrelation  [II]  suggests  the  existence 
of  a  spatial  counterpart  to  dispersion  cancellation.  How¬ 
ever,  no  experimental  observation  of  a  spatial  dispersion- 
cancellation  effect  has  been  reported  so  far. 

The  longitudinal  component  of  the  wave  vector,  on  the 
other  hand,  sets  up  the  phase-matching  relation  that  estab¬ 
lishes  conditions  for  an  effective  energy  conversion  be¬ 
tween  three  interacting  waves,  pump,  signal,  and  idler. 
Since  the  longitudinal  wave  vector  depends  both  on  fre¬ 
quency  and  on  transverse  momentum,  this  condition  sets  a 
specific  relation  between  the  frequency  and  the  emission 
angle  of  down-converted  radiation.  In  other  words,  the 
frequency  and  transverse  momentum  degrees  of  freedom 
cannot  be  factorized  and  the  overall  quantum  state  is  con¬ 
currently  entangled  in  both  o)  and  q  (multiparameter  en¬ 
tanglement).  This  leads  to  several  interesting  effects  where 
the  manipulation  of  a  spatial  variable  affects  the  shape  of 
the  temporal  interference  pattern  and  also  polarization 
interference  pattern  [12]  . 

In  this  Letter,  we  exploit  the  multiparameter  en¬ 
tangled  states  generated  by  SPDC  to  demonstrate  the 
effect  of  even-order  spatial  aberration  cancellation.  We 
use  an  SPDC  source  to  produce  momentum-anticorrelated 
photons  and  we  modulate  their  wave  fronts  by  a  transfer 
function  H(q).  Because  of  the  correlations  between  q  and 
0),  the  manipulation  in  the  q  space  will  introduce 
changes  in  the  o)  space.  Therefore  the  spatial  wave  front 
modulation  will  affect  the  temporal  interference  pat¬ 
tern,  which  can  be  observed  using  a  polarization  two- 
photon  interferometer  [12,13].  With  this  technique  we 
show  that,  due  to  the  anticorrelation  of  the  transverse 
momenta,  the  even-order  aberrations  are  canceled  out, 
and  only  the  odd-order  contributions  influence  the 
resulting  interference  pattern.  We  believe  this  effect  may 
lead  to  interesting  applications  in  the  field  of  quantum 
imaging. 
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The  setup  of  our  experiment  is  sketched  in  Eig.  1 .  A  laser 
diode  with  a  single  longitudinal-mode  selection  (405  nm, 
35  mW)  pumps  a  L5-mm  thick  BBO  (BaB204)  crystal  that 
is  cut  for  a  collinear  degenerate  type-II  phase  matching. 
Approximating  the  pump  beam  with  a  plane  wave,  the  two 
orthogonally  polarized  photons  emitted  by  the  crystal  can 
be  described  by  the  quantum  state  [10]: 

\ip)  =  j dq  j dco^iq,  (o)d\{q,  Hq  +  (o) 

X  (-q,  Ho  -  «)|0),  (1) 


two  distinct  paths,  one  towards  a  flat  mirror  (EM)  and  the 
other  towards  a  deformable  mirror  (DM).  Each  photon 
passes  through  a  4-/  system  comprising  a  lens  (LI)  of 
focal  length  /  =  200  mm  positioned  at  a  distance  /  from 
the  output  plane  of  the  crystal,  and  the  same  distance  / 
from  the  mirror.  On  the  way  from  the  crystal  to  the  mirror, 
the  lens  maps  each  wave  vector  component  to  a  different 
point  x(q,  co)  on  the  mirror  surface.  The  limited  down- 
conversion  bandwidth  (about  30  nm  for  a  collection  angle 
of  25  mrad)  allows  us  to  neglect  the  frequency  dependence 
of  the  4-/  system;  assuming  that  the  lens  is  achromatic: 


where 


^(q,  (o)  =  sinc^^^^^'  j^4^A(q,ftj)/2)  ^2) 

The  phase-mismatch  function  A(q,  co)  in  the  paraxial 
approximation  takes  the  form 

A(q,  co)  =  —  coD  +  M^2  '  Q  ,  (3) 

kp 

where  D  is  the  difference  between  the  inverse  of  the  group 
velocities  of  the  ordinary  and  the  extraordinary  waves 
inside  the  birefringent  crystal  and  M  is  their  spatial  walk- 
off  in  the  vertical  direction  62-  The  term  is  due  to 

diffraction  during  the  propagation  of  photons  through  the 
crystal.  In  case  of  a  BBO  crystal,  phase-matched  for  a 
degenerate  (Aq  =  810  nm)  collinear  type-II  down- 
conversion,  D  =  182  ps/mm  and  M  =  0.0723. 

A  polarizing  beam  splitter  (PBS)  separates  the  horizon¬ 
tally  polarized  photon  and  the  vertically  polarized  one  into 


FIG.  1  (color  online).  Schematic  of  the  experimental  setup. 
Examples  of  aberrations  induced  by  the  deformable  mirror  are 
illustrated  in  the  inset:  the  even-parity  aberrations  on  the  left  are 
canceled,  while  the  odd-parity  ones  on  the  right  alter  the  shape  of 
the  interference  pattern. 


x(q. ")  =  -rAq  (4) 

k{0))  kg 

The  deformation  of  the  mirror  surface  at  point  x  can  be 
described  by  the  function  ^(x).  A  photon  focused  to  the 
point  X  by  the  lens  will  travel  a  distance  ^(x)  to  the  mirror 
surface,  will  be  reflected  and  will  travel  a  distance  ^(x) 
back  to  the  lens  focal  plane.  Therefore  it  will  acquire  a 
phase  shift: 

cp{x)  -  2kM,  (5) 

After  reflection  from  the  mirror,  the  same  lens  maps  every 
point  back  to  a  wave  vector.  Mathematically,  the  trans¬ 
formation  induced  by  the  4-/  system  can  be  described  by 
the  transfer  function: 

^(q)  =  (6) 

where  the  pupil  function  p{x)  describes  the  circular  aper¬ 
ture  of  the  mirror. 

The  deformable  mirror  [14]  consists  of  a  thin  nitrocel¬ 
lulose  silver-coated  membrane  (12  mm  diameter,  5  /mm 
thick,  initial  flatness  less  than  20  nm  rms)  that  is  deformed 
by  electrostatic  forces  created  when  a  voltage  drop  (maxi¬ 
mum  270  V)  is  applied  to  37  electrodes.  The  action  pro¬ 
duced  by  each  actuator  was  mapped  by  measuring  the 
induced  deformation  with  a  Zygo  interferometer,  creating 
an  influence  function  matrix. 

In  addition,  each  photon,  traveling  from  the  PBS  to  the 
mirror  and  back  passes  twice  through  a  quarter- wave-plate 
oriented  at  45°  and  flips  its  polarization  state.  This  way  the 
photon  that  has  been  transmitted  is  now  reflected  at  the 
polarizing  beam  splitter  (PBS)  and  vice  versa,  resulting  in 
both  photons  leaving  the  modulation  section  together  to¬ 
wards  the  polarization  interferometer. 

The  polarization  interferometer  [13,15]  consists  of  a 
birefringent  delay  line  (DL,  made  of  two  sliding  quartz 
wedges)  providing  a  variable  delay  r  and  a  nonpolarizing 
beam  splitter  (BS)  that  splits  the  photons  in  two  paths 
directed  to  two  single-photon  detectors  Di  and  £>2.  A 
polarizer  oriented  at  45°  is  placed  in  front  of  each  detector 
in  order  to  erase  information  about  the  polarization  of  the 
incoming  photon.  Photons  were  collected  by  a  lens  in  each 
arm  and  focused  onto  the  detector’s  active  area.  To  max¬ 
imize  the  spatial  collection  capability  we  used  two  open- 
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face  (180  yam  diameter)  single-photon  silicon  avalanche 
photodiodes.  Using  a  fiber  coupler  would  limit  the  number 
of  collected  spatial  modes.  All  experiments  have  been 
performed  with  an  8-mm  diameter  pinhole  placed  at 
330  mm  from  the  output  plane  of  the  4-/  system,  therefore 
collecting  photons  from  an  angle  of  about  25  mrad.  A 
dichroic  mirror  and  a  pair  of  interference  bandpass  filters 
with  a  bandwidth  that  is  greater  than  that  of  down- 
conversion  have  been  used  to  reject  the  residual  pump 
radiation  and  the  background  light.  The  number  of  coinci¬ 
dences  acquired  as  a  function  of  optical  polarization  delay 
T  shows  a  high- visibility  quantum  interference  pattern  in 
the  form  of  a  dip  [12]. 

Since  the  photon-counting  detectors  are  slow,  compared 
with  the  coherence  time  of  down-converted  photons,  and 
since  their  surface  is  larger  than  the  spot  size,  the  expres¬ 
sion  for  the  coincidence  rate,  in  the  paraxial  approxima¬ 
tion,  is  [16] 

Rc(t)  =  ^o[i-a(i-;^)wm(t)],  (7) 

where  Rq  is  the  background  coincidence  rate,  A(Qf)  is  a 
triangular  function  [A(af)  =  1  —  |Qf|  if  |Qf|<l,  and 
A(6r)  =  0,  otherwise],  and 

+  q']sinc|^MLe2  •  (q  +  qOA^l  “  j 

X  (g) 

The  function  A(1  —  represents  a  usual  triangular  dip 
one  obtains  in  type-II  quantum  interferometry  when  work¬ 
ing  in  the  single  spatial-mode  approximation  (using  narrow 
pinholes).  The  function  Wm{t)  takes  into  account  the 
deformation  of  the  triangular  dip  induced  by  the  modula¬ 
tion  of  the  transverse  wave  vectors  and  the  Fourier  trans¬ 
form  of  the  shape  of  the  detection  apertures  PaLq]* 
particular,  the  function  Wm{t)  describes  how  manipulation 
in  the  q  space  by  a  filter  //(q)  is  converted  into  a  modifi¬ 
cation  of  the  temporal  interference  pattern,  by  means  of  the 
coupling  between  wave  vectors  and  frequencies  set  by  the 
phase-matching  conditions. 

If  the  detection  apertures  are  sufficiently  large,  the  func¬ 
tion  Pa[q  +  can  be  well  approximated  by  a 
delta  function,  so  that  Eq.  (8)  can  be  simplified  to 


If  the  function  //(q)  has  a  circular  symmetry,  its  phase 
^(q)  =  arg{//(q)}  can  be  expanded  using  a  set  of  Zernike 
polynomials,  which  are  orthogonal  on  the  unit  circle  [17], 

•PCq)  =  Z  S  VnmRnip)  COs{md),  (10) 


where  q  =  (p  cos^,  p  sin^)  and  m  =  —n,  —n  +  2,  —n  + 
A, ..  .,n.  Using  the  fact  that  —  q  =  [p cos(^  +  tt), 
p  sin(^  +  tt)],  and  that  if  m  is  even,  then  cos[m(^  +  tt)]  = 
cos(m^),  while  if  m  is  odd,  then  cos[m(^  +  tt)]  = 
—  cos(m^),  one  gets 

^(q)  -  ^(-q)  =  2^  X  (PnmRnip)cos{md).  (11) 

n  m  odd 

This  means  that  only  Zernike  polynomials  with  m  odd 
(and  consequently  n  odd)  will  contribute  to  the  shape  of  the 
interference  pattern.  For  example,  contributions  from  as¬ 
tigmatism  (n  =  2,  m  =  ±2),  defocus  (n  =  2,  m  =  0),  and 
spherical  aberration  (n  =  4,  m  =  0)  will  all  be  canceled. 
On  the  contrary,  coma  (n  =  3,  m  =  ±1)  and  trefoil  (n  = 
3,  m  =  ±3)  will  be  present. 

We  studied  different  types  of  aberrations:  coma  along 
the  V  axis  (n  =  3,  m  =  + 1),  coma  along  the  y  axis  (n  =  3, 
m  =  —  1),  astigmatism  (n  =  2,  m  =  2),  and  the  aberration 
corresponding  io  n  =  A  and  m  =  A.  The  experimental  data 
for  coma  oriented  along  the  v  and  y  (parallel  to  the  vertical- 
polarization)  directions  are  presented  in  Fig.  2.  Because  of 
the  multiparameter  entanglement  of  the  two-photon  state, 
the  wave  front  distortion  induced  by  the  deformable  mirror 
modulates  the  spectral  degree  of  freedom,  resulting  in  a 
modification  of  the  temporal  interference  pattern.  We  in¬ 
creased  the  maximum  amplitude  of  the  mirror  deformation 
from  0.2  yam  (—0.25 Aq)  to  0.75  yam  (~  Aq)  (peak-to  valley 
deformation)  for  coma  along  v  and  from  0.2  yam  to 
0.4  yam  (— 0.5Ao)  for  coma  along  y  axis.  The  shape  of 
the  interference  pattern  is  changed  dramatically  when  the 
intensity  of  the  coma  aberration  increases.  Theoretical 
predictions  (solid  lines)  based  on  the  experimental  parame- 


FIG.  2  (color  online).  Coincidence-rate  interference  patterns 
when  coma  (n  =  3,m=±l)  along  the  x  axis  (on  the  left)  and 
along  the  y  axis  (on  the  right)  is  imposed  on  the  deformable 
mirror.  Solid  lines  illustrate  theoretical  fitting  with  experimental 
parameters.  The  initial  relative  tilt  between  the  deformable 
mirror  (DM)  and  the  flat  mirror  (FM)  is  used  as  an  adjustable 
parameter  to  account  for  the  imperfectness  of  experimental 
alignment  between  two  arms.  The  shapes  of  adaptive  mirror 
deformation  producing  selected  aberrations  are  illustrated  in  the 
insets. 


233603-3 


PRL  101,  233603  (2008) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
5  DECEMBER  2008 


1.2 


0.8 


0.6 


^  0.4 


deloy  (fs)  delay  (fs) 


FIG.  3  (color  online).  Interference  patterns  in  the  coincidence 
rate  when  astigmatism  at  45°  (left  picture)  and  the  aberration 
described  by  n  =  A,  m  =  A  are  dialed  on  the  mirror.  The 
intensity  of  the  aberrations  has  been  scanned  from  0.2  to 
0.8  yum  peak-to-valley  deformation.  The  interference  pattern  is 
insensitive  to  the  aberrations  within  experimental  errors. 


ters  are  superposed  to  experimental  data  in  Eig.  2:  the 
matching  between  the  two  curves  is  pretty  good. 

Experimentally  obtained  data  for  astigmatism  (with 
symmetry  axes  oriented  at  45°  with  respect  to  the 
x-y  axes)  and  for  the  aberration  identified  hy  n  =  A  and 
m  =  A  are  shown  in  Eig.  3.  The  effect  of  these  even-order 
aberrations  is  effectively  canceled  out  due  to  the  spatial 
correlations  between  the  photons  in  parametric  down- 
conversion.  Therefore,  such  type  of  spatial  aberrations  do 
not  affect  the  shape  of  the  dip  so  that  the  known  quantum 
interference  pattern  is  retained. 

This  effect  has  a  clear  analogy  with  even-order  fre¬ 
quency  dispersion  cancellation  due  to  frequency  entangle¬ 
ment  in  SPDC.  In  case  of  spectral  dispersion,  the  use  of  a 
nonmonochromatic  pump  reduces  the  degree  of  correlation 
between  spectral  components  of  entangled  photons  and 
degrades  the  dispersion  cancellation  effect.  In  our  case  of 
even-order  aberration  cancellation,  sharp  spatial  correla¬ 
tions  can  be  obtained  only  in  the  approximation  of  a  plane- 
wave  pump  beam.  Therefore,  wave  vector  correlations  get 
weaker  for  focused  pump  beams  and  the  aberration  can¬ 
cellation  effect  also  degrades  when  the  pump  beam  is 
focused  tightly  on  the  crystal.  Eurthermore,  the  spectral 
dispersion  cancellation  effect  works  in  the  limit  of  slow 
detectors  because  it  requires  integration  over  the  temporal 
degree  of  freedom.  Similarly,  the  even-order  aberration 
cancellation  works  well  only  in  the  situation  when  the 
collection  apertures  used  in  experiment  are  sufficiently 
large  to  enable  effective  integration  over  the  spatial  degrees 
of  freedom  [16]. 

The  question  remains  whether  the  effect  we  have  re¬ 
ported  is  purely  quantum,  or  some  classical  counterpart  can 
be  envisioned,  as  in  the  case  of  spectral  dispersion  cancel¬ 
lation.  We  believe  that  a  classical  optical  configuration 
mimicking  even-order  aberration  cancellation  could  poten¬ 


tially  be  discovered  by  exploiting  a  light  source  with  strong 
degree  of  spatial  intensity  correlation  similar  to  optical 
speckles. 

In  conclusion,  we  have  demonstrated  experimentally  the 
effect  of  cancelling  even-order  spatial  aberration  in  quan¬ 
tum  interferometry  using  entangled-photon  states  gener¬ 
ated  in  a  type-II  spontaneous  parametric  down-conversion 
process.  This  effect  may  prove  helpful  in  enhancing  the 
spatial  resolution  in  quantum  imaging. 
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verse  wave  vectors  in  type-1 1  parametric  down-conversion  using  adaptive  optics  techniques.  In  particular,  we 
discuss  the  possibility  of  spatial  walk-off  compensation  in  quantum  interferometry  and  an  effect  of  even-order 
spatial  aberration  cancellation. 

DOI:  10.1103/PhysRevA. 79.062304  PACS  number(s):  03.67.Bg,  42.50.St,  42.50.Dv,  42.30.Kq 


I.  INTRODUCTION 

Quantum  entangiement  [1]  is  a  valuable  resource  in  many 
areas  of  quantum  optics  and  quantum  information  process¬ 
ing.  One  of  the  most  widespread  techniques  for  generating 
entangied  optical  states  is  spontaneous  parametric  down- 
conversion  (SPDC)  [2-5].  SPDC  is  a  second-order  noniinear 
optical  process  in  which  a  pump  photon  is  spiit  into  a  pair  of 
new  photons  with  conservation  of  energy  and  momentum. 
The  phase-matching  reiation  estabiishes  conditions  to  have 
efficient  energy  conversion  between  the  pump  and  the  down- 
converted  waves,  caiied  signal  and  idler.  This  condition  also 
sets  a  specific  reiation  between  the  frequency  and  the  emis¬ 
sion  angie  of  down-converted  radiation.  In  other  words,  the 
quantum  state  emitted  in  the  SPDC  process  cannot  be  factor¬ 
ized  into  separate  frequency  and  wave-vector  components. 
This  leads  to  several  interesting  effects  where  the  manipula¬ 
tion  of  a  spatial  variable  affects  the  shape  of  the  polarization- 
temporal  interference  pattern.  For  example,  the  dependence 
of  polarization-temporal  interference  on  the  selection  of  col¬ 
lected  wave  vectors  was  studied  in  detail  in  [6]. 

Here  we  engineer  the  quantum  state  in  the  space  of  trans¬ 
verse  momentum  and  we  study  how  this  spatial  modulation 
is  transferred  to  the  polarization-spectral  domain  by  means 
of  quantum  interferometry.  We  will  focus  on  type-1 1  SPDC 
using  birefringent  phase  matching  since  the  correlations  be¬ 
tween  wave  vectors  and  spectrum  are  stronger  than  those 
employing  other  phase-matching  conditions. 

Our  aim  is  twofold.  From  one  point  of  view,  we  study  the 
effect  of  spatial  modulations  on  temporal  quantum  interfer¬ 
ence.  This  could  be  useful,  for  example,  in  quantum  optical 
coherence  tomography  (QOCT)  [7,8].  When  focusing  light 
on  a  sample  with  nonplanar  surface,  the  photons  will  acquire 
a  spatial  phase  distribution  in  the  far  field,  which  may  per¬ 
turb  the  shape  of  the  interference  dip.  Our  results  will  pro¬ 
vide  a  tool  to  understand  this  effect. 

From  a  second  point  of  view,  we  would  like  to  study  and 
characterize  spatial  modulation  as  a  tool  for  quantum  state 
engineering.  This  may  find  application  in  the  field  of  quan¬ 
tum  information  processing,  where  it  is  important  to  gain  a 
high  degree  of  control  over  the  production  of  quantum  en¬ 
tangled  states  entangled  in  one  or  more  degrees  of  freedom 
(hyperentanglement). 

We  start  by  introducing  a  theoretical  model  of  a  type-1 1 
quantum  interferometer,  comprising  the  polarization,  spec¬ 


tral,  and  spatial  degrees  of  freedom  (Sec.  II).  A  modulation  in 
the  wave-vector  space  is  provided  by  an  adaptive  optical 
setup  and  equations  for  the  polarization-temporal  interfer¬ 
ence  pattern  in  the  coincidence  rate  are  derived.  In  Sec.  Ill, 
we  introduce  a  numerical  approach  for  practical  evaluation 
of  the  results  of  the  theoretical  model,  discussing  a  few  ex¬ 
amples  for  general  spatial  aberrations. 

In  Sec.  IV  we  will  highlight  and  discuss  theoretically  two 
interesting  special  cases.  The  first  one  is  the  possibility  of 
restoring  high  visibility  in  type-1 1  quantum  interference  with 
large  collection  apertures.  In  some  situations,  to  collect  a 
higher  photon  flux  or  a  broader  photon  bandwidth,  it  can  be 
useful  to  enlarge  the  collection  apertures  of  the  optical  sys¬ 
tem.  But  when  dealing  with  type-1 1  SPDC  in  birefringent 
crystals,  for  large  collection  apertures  the  effect  of  spatial 
walkoff  introduces  distinguishability  between  the  photons, 
leading  to  a  reduced  visibility  of  temporal  and  polarization 
quantum  interference.  We  will  show  that  high  visibility  can 
be  restored  with  a  linear  phase  shift  along  the  vertical  axis. 

The  second  effect  is  the  spatial  counterpart  of  spectral 
dispersion  cancellation  [9,10].  In  the  limit  of  large  detection 
apertures,  the  correlations  between  the  photons'  momenta 
will  cancel  out  the  effects  of  even-order  aberrations,  exactly 
as  in  the  limit  of  slow  detectors  the  frequency  correlations 
cancel  out  the  even-order  terms  of  spectral  dispersion.  The 
experimental  demonstration  of  this  effect  has  been  reported 
recently  [11]. 

As  we  proceed  from  the  general  case  of  Secs.  II  and  III 
into  the  specific  examples  of  Sec.  IV,  we  will  gradually  see 
that  optical  aberration  is  a  subject  with  two  very  different 
faces.  On  one  hand,  aberrations  in  optical  components  are 
normally  seen  as  undesirable,  si  nee  they  lead  to  distortions  in 
imaging.  We  will  see  that  these  unwanted  spatial  modula¬ 
tions  may  to  some  extent  be  canceled.  On  the  other  hand,  we 
will  find  that  such  spatial  modulations  may  also  be  turned 
into  a  useful  tool:  by  deliberately  introducing  spatial  modu¬ 
lations  (in  effect,  purposely  adding  aberrations  in  a  con¬ 
trolled  manner),  we  can  produce  useful  effects  such  as  the 
restoration  of  visibility  mentioned  above.  The  examples  we 
provide  in  Sec.  IV  will  illustrate  these  two  aspects  and  will 
show  that  both  can  benefit  from  more  detailed  study  of  the 
interplay  between  spatial  modulations  and  spectral  correla¬ 
tions. 
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FIG.  1.  (Color  online)  Scheme  of  the  proposed  setup.  Horizon¬ 
tally  polarized  photons  from  type-1 1  SPDC  are  assigned  a  phase 
dependent  on  the  photon  transverse  momentum  ^c>(Qo)-  ''''hUe  the 
vertically  polarized  ones  are  assigned  a  phase  <^e(Qe)-  The  modu¬ 
lated  photons  enter  a  type-1 1  quantum  interferometer,  which  records 
the  coincidence  count  rate  as  a  function  of  the  delay  t  between  the 
photons  given  by  an  appropriate  delay  line. 


II. THEORETICAL  MODEL 

Consider  the  scheme  in  Fig.  1.  A  laser  beam  pumps  a 
nonlinear  material  phase  matched  for  type-1 1  parametric 
down-conversion,  creating  a  pair  of  entangled  photons.  Each 
of  the  generated  photons  passes  through  a  Fourier-transform 
system,  and  then  enters  a  modulation  system  which  trans¬ 
forms  transverse  wave  vectors  for  the  horizontally  (H)  polar¬ 
ized  photon  according  to  the  transfer  function  G  i(qi)  and  for 
the  vertically  (V)  polarized  photon  according  to  G2(q2).  Af¬ 
ter  being  modulated  In  q  space,  the  photons  enter  a  type-1 1 
Interferometer.  A  nonpolarizing  beam  splitter  (BS)  creates 
polarization  entanglement  from  the  polarization-correlated 
pair  emitted  by  the  source.  The  beams  at  the  output  ports  of 
the  beam  splitter  are  directed  toward  two  single-photon  de¬ 
tectors.  Two  polarizers  at  45°  before  the  BS  restore  Indlstin- 
guishability  in  the  polarization  degree  of  freedom.  An  adjust¬ 
able  delay  line  r  is  scanned  and  the  coincidence  rate  R(t) 
between  the  detection  events  of  the  two  detectors  is  recorded. 
An  aperture  is  placed  before  the  beam  splitter  to  select  an 
appropriate  collection  angle. 


A.  Notation 


Consider  a  monochromatic  plane  wave  of  complex  ampli¬ 
tude  E  (r)  =£  with  r=(x,y,z).  For  a  given  wavelength 
corresponding  to  a  frequency  fl,  the  wave  vector  can  be 
split  into  a  transverse  component  q=(gx-qy)  and  a  longitudi¬ 
nal  component  /3(q;fl): 

k=[q,yS(q;a)].  (1) 


The  wave  number  Is 


fc(a)  = 


c 


(2) 


The  longitudinal  component  of  the  wave  vector  is 

Mq-il)  =  Mll)-|qp.  (3) 

Therefore  the  electric  field  at  the  position  r  and  time  f  can  be 
written  as 


f(r;t) 


dn£(q,n)e-'‘'-''e-''^(‘''^V^f, 


(4) 


where  p=(x,y). 


In  paraxial  approximation  |q|^<^/c^(fl),  so  that 

For  a  quasimonochromatic  wave  packet  centered  around  the 
frequency  Ho,  one  can  write  il=dQ+v,  with  This 

expression  can  be  approximated  by 

yS(q,a)-ko  +  ;f-^-  (6) 

Uq  ZKq 

where  kQ=k(ilo)  and  Uo=(  is  the  group  velocity 

for  the  propagation  of  the  wave  packet  through  the  material. 

B.  State  generation 

Using  first-order  time-dependent  perturbation  theory,  the 
two-photon  state  at  the  output  of  the  nonlinear  crystal  can  be 
calculated  as 


ItA) (t)|0),  (7) 

where  the  interaction  Hamiltonian  is 

H,(t)  =  ^J  (8) 

The  strong,  undepleted  pump  beam  can  be  treated  classically. 
Assuming  a  monochromatic  plane  wave  propagating  along 
the  z  direction. 


Ep(r,t)=Epe'TV-V).  (9) 

The  signal  and  Idler  photons  are  described  by  the  following 
quantum  field  operators: 

Ej"V-t)=  j  dqjj  dwje4A(‘ii'"i>^+‘ii'P-"ida(q^.,w^.),  (lO) 
where  j=e,o. 

The  biphoton  quantum  state  at  the  output  plane  of  the 
nonlinear  crystal  is  [12] 


di^(q,  v)ao(q,no  +  v)al(-q,flo  -  v)|0). 


(11) 


Two  photons  are  emitted  from  the  nonlinear  crystal,  one 
horizontally  polarized  (ordinary  photon)  and  the  other  verti¬ 
cally  polarized  (extraordinary  photon),  with  anti  correlated 
frequencies  and  emission  directions. 

In  the  case  of  a  single  bulk  crystal  of  thickness  L  and 
constant  nonlinearity  Xo,  the  probability  amplitude  for  having 
the  signal  photon  in  the  mode  (q,flo+i')  and  the  idler  in  the 
mode  (-q,fio-v)  is 


<l>(q,i')  =sinc 


tA(q,v) 

2 


g/[A(q,>')L/2] 


(12) 


For  type-1 1  col  I  inear  degenerate  phase  matching,  the  phase- 
mismatch  function  A(q,  v)  can  be  approximated  to  be 
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A(q,v)=-vD -I-Me2'q  +  -— ,  (13) 

Kp 

where  D  is  the  difference  between  the  inverse  of  the  group 
veiocities  of  the  ordinary  and  extraordinary  photons  inside 
the  birefringent  crystai  and  the  quadratic  term  in  q  is  due  to 
diffraction  in  paraxiai  approximation.  The  remaining  term  is 
the  first-order  approximation  for  the  spatiai  waikoff. 

C .  Propagation 

Consider  a  photon  described  by  the  operator  aj(q,fl)  (po- 
iarization  j  =e,o,  frequency  fl,  and  transverse  momentum  q). 
Its  propagation  through  an  optical  system  to  a  point  on  the 
output  piane  is  described  by  the  optical  transfer  function 
H j(Xi(,q;fl).  In  our  setup,  the  field  at  the  detector  will  be  a 
superposition  of  contributions  from  the  ordinary  and  extraor¬ 
dinary  photons.  The  quantized  electric  fields  at  the  detector 
planes  are 

^A'^\xA,tA)=J  dqj  dcoe'‘^^^[He(XA,q;(o)ae(q,co) 

+  Ho(XA,q;(o)aJq,(o)], 

E^B^\xB,tB)=J  dqj  d«e'“'e[He(XB,q;w)ae(q,«) 

-l-Ho(Xs,q;w)ao(q,«)].  (14) 

The  probability  amplitude  for  detecting  a  photon  pair  at  the 
detector  planes,  with  space-time  coordinates  (XA,tA)  and 
(Xe-te),  is 


A(XA,XB;tA,tB)  =(0|£Y^(x,,t,)E«(XB,ts)|*A).  (15) 

For  the  bi photon  probability  amplitude  we  get 

A{XA,XB;tA,tB)=  j  dqodqed(Oodcog(l>(qo,qg;(Oo,cog) 

X  [H  e(XA  -  qe  I  o(Xb  ,  qo  I 

+  H  o(X;, ,  qo ;  Wo)H  e(XB ,  qe;  . 

(16) 

This  probability  amplitude  represents  the  superposition  of 
two  possible  events  leading  to  a  coincidence  count  in  the 
detectors: 

(1)  the  V  polarized  photon  with  momentum  qg  and  fre¬ 
quency  cog  going  through  the  lower  branch  to  arrive  at  posi¬ 
tion  X;,  in  detector  A,  while  the  H  polarized  photon  with 
momentum  q^  and  frequency  wq  goes  through  the  upper 
branch  to  arrive  at  position  Xb  in  detector  B;  and 

(2)  the  V  polarized  photon  with  momentum  qg  and  fre¬ 
quency  cog  going  through  the  lower  branch  to  arrive  at  posi¬ 
tion  Xb  in  detector  6,  while  the  H  polarized  photon  with 
momentum  q^  and  frequency  wg  goes  through  the  upper 
branch  to  arrive  at  position  Xa  in  detector  A. 

Since  the  superposition  is  coherent,  there  are  quantum- 
interference  effects  between  the  two  probability  amplitudes. 


1.  State  engineering  section 

In  the  state  engineering  section,  each  of  the  two  branches 
consists  of  a  pair  of  achromatic  Fourier-transform  systems 
coupled  by  a  spatial  light  modulator  or  a  deformable  mirror. 
Each  Fourier-transform  system  consists  of  a  single  lens  of 
focal  length  f,  separated  from  the  optical  elements  before 
and  after  it  by  a  distance  f.  The  first  Fourier  system  maps 
each  incident  transverse  wave  vector  q  on  the  plane  Il/n  to  a 
point  x(q)  on  the  Fourier  plane  11/:: 

x(q)  =  7^q,  ko  =  ^,  (17) 

Kq  c 

where  f  is  the  focal  length  of  the  Fourier-transform  system. 
Since  we  assume  that  the  system  is  achromatic  for  a  certain 
bandwidth  around  a  central  frequency  Hq,  the  position  x(q) 
depends  only  on  q  and  not  on  w. 

The  spatial  modulator  assigns  a  different  amplitude  and 
phase  to  the  light  incident  on  each  point,  as  described  by  the 
function  G (x)=t(x)e'^\  Each  point  is  then  mapped  back  to 
a  wave  vector  on  the  plane  rtgut  by  the  second  achromatic 
Fourier-transform  system. 

Using  the  formalism  of  Fourier  optics  [13],  the  transfer 
function  between  the  planes  !!/„  and  Ilgyt  can  be  calculated 
to  be 


hi(Xi,X3)=  j  dx6(x)e-'(*^o/0x-(xi+x3)^ 

The  corresponding  momentum-transfer  function  is 

f 


Hi(qi,q3)=G 


-qi 


^qi-qs)- 


(18) 


(19) 


2.  interferometer 

After  the  plane  Flout  the  two  photons  enter  a  type-1 1  quan¬ 
tum  interferometer.  Each  propagates  in  free  space  to  a  bire¬ 
fringent  delay  line  and  a  detection  aperture  p(x)  to  be  finally 
focused  to  the  detection  planes  by  means  of  lenses  of  focal 
length  fg.  Following  the  derivation  in  [6]  the  transfer  func¬ 
tion  is 


H2(Xi,q;w) 


I 


h(xi,Xj;(o)e'^dXi 


_  g/(w/C)(cfi+Cf2+/'o)  0xp 


.‘^IXiP 


^g-/(ccli/2a))|qpp 


(20) 


where  P(q)  is  the  Fourier  transform  of  |p(x)p. 

A  combination  of  the  two  different  stages  is  described  by 
the  transfer  function 


H2(Xj,q„;«), 


(21) 


where  the  two  functions  Gi(q)  and  G2(q)  are  the 
momentum-transfer  functions  which  describe  the  modulation 
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imparted,  respectively,  on  the  ordinary  and  the  extraordinary 
photons. 


D.  Detection 

Since  the  singie-photon  detectors  used  in  quantum  optics 
experiments  are  siow  with  respect  to  the  temporal  coherence 
of  the  photons  and  their  area  is  iarger  than  the  spot  onto 
which  the  photons  are  focused  by  the  coilection  lens,  we 
integrate  over  the  spatial  and  temporal  coordinates.  There¬ 
fore  the  coincidence  count  rate  expressed  in  terms  of  the 
biphoton  probability  amplitude  is 

R(t)=  j  dXA  j  dxsj  dtA  j  dtBlA(XA,XB;tA,tB  f.  (22) 

Following  the  derivation  described  in  Appendix  A,  one 
gets 


R(t)=Ro 


1-AI 


1 


where  A(x)  is  the  triangular  function 


A(x)  = 


|x|<l 

|X|>1. 


(23) 


(24) 


Therefore,  the  coincidence  count  rate  fifr)  is  given  by  the 
summation  of  a  background  level  Rq  and  an  interference  pat¬ 
tern  given  by  the  trianguiar  dip  A(l-|^)  that  one  gets  when 
working  with  narrow  apertures,  moduiated  by  the  function 
Wm(t)  which  depends  on  the  detaiis  of  the  adaptive  optical 
system. 

The  expressions  for  Rq  and  WM(r)  are 


FIG.  2.  (Color  online)  Example  of  the  numerical  approach 
adopted  to  evaluate  Eqs.  (25)  and  (26).  The  spatial  modulation  sur¬ 
face  is  discretized  in  sufficiently  small  squares  over  which  the  phase 
is  averaged. 

III.  NUMERICAL  SOLUTIONS  FOR  A  GENERAL 
PHASE  SHIFT 


Ro  = 


dq'  sinc[MLe2'(q-q')] 


xG 


xG 


g-/(M  L  /2)e2-(q-q  ')e/(2c/i/'fp)[|qP-|q '  P] 


xPA(q-q')PB(-q+q') 


(25) 


and 


dq'  sine 


MLej-  (q  -Fq')A 


ko  / 


xP,[q+q']PB[-(q+q')].  (26) 

In  the  following  we  will  assume  there  is  spatial  modula¬ 
tion  only  on  one  of  the  photons;  therefore  we  set  G2(q)  =  l. 


Numerically  solving  for  the  quantities  in  Eqs.  (25)  and 
(26)  in  the  case  of  a  general  aberration  may  be  computation¬ 
ally  demanding.  Here,  we  propose  an  approximation,  valid  in 
the  case  where  the  function  G  (x)  changes  smoothly  over  the 
mirror  surface,  as  it  is  the  case  in  experimentally  relevant 
situations.  This  model  is  also  interesting  from  the  practical 
point  of  view,  since  in  many  cases  adaptive  optical  systems 
are  implemented  using  spatial  light  modulators  or  segmented 
deformable  mirrors,  where  the  modulation  surface  is  parti¬ 
tioned  into  small  pixels. 

Suppose  we  partition  the  Fourier  plane  IT/:  into  small 
squares  (pixels)  of  side  d  (Fig.  2).  Let  us  define  the  rectan¬ 
gular  function 


n(x)  = 


0  if  |x|  >  I 
1  if  |x|<|. 


(27) 


The  pixel  (l,m)  is  identified  by 


ai,rr,(x,y)=n 


X  , 

n 

y 

+ 

^  +m 

d 

d 

(28) 


selecting  the  area 
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(m-l)d<y<(m  +  l)d.  (29) 

We  approximate  the  value  of  the  phase  in  each  square  by 
the  mean  value  of  the  actual  phase  within  the  square: 


■ii 


<p;m  =  :^  dxdy<p(x,y)n 


n 


y 

^  +m 
d 


That  is, 


gi<p(x,y)  ^  ^  gi<pij„n[(x/d)+l]n[(y/d)+m]^ 
l,m 


In  this  case  (see  Appendix  B  for  a  justification) 


(31) 
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^  g/<A(,mn[(x/c/)+/]n[(y/c/)+m]  _  ^  e'<A(,mn 
l,m  l,m 


X  , 

n 

y 

+ 

^  +m 

d 

d 

(32) 


(30)  Substituting  this  expression  into  Eq.  (23)  and  collecting  the 
integrations,  one  finds 


R (r)  =  E  E  e-'>'--^^.Ma„/„^(r),  (33) 

l,m  X,fJL 


where 


«;x 


=  j  dqj  dO,n 


n 


Qii2dykp)iql-Ql)P[q,+Q,V[-(q,+Q,)] 


and 


Irn.fiiT) 


=  j  dQyj  dOyU 


Id^y-^ 


n 


gj(2di/kp}(qj-Qpg-j(M/DHqy-0y) 


xSinc 

Performing  the  integrations  one  gets 

aa  =  J  dxP(x)A 
and 


ML(g.-EQ.)A  1- 


2t 


P[gy  +0y]P[-(gy  +Qy)]. 


a-'--' 


•  A 

smcf  —^xA 


q/(c/c/i/0(/-x)x 


-\ 


lm^(r)=  I  dxP(x)A 


fx 

kd 


-{m+fi) 


sinci  M  Lx  A 


1- 


2t 


.  \2kd(2di 

inr<  -  - y  — 


xsinci 


M 


I  f 


fx  , 


g/(M/f)[(2c/i/J:p)x-(M  /D)T\(m-^0  _ 


A  similar  expression  can  be  found  for  the  background  coincidence  rate: 


l,m  k,jui 


where 


and 


-! 

R(yl  =  l  dxP(x)A 


rW=  dxP(x)A 


.  I 2ddi 
sincf  —^xA 


f  2ddi 

sinc(MLx)sincj  —^xA 


fx 


q/(c/c/i/0(/+x)x 


Q/[(cfcfi/0(m+/x)-M/./2]x 


(34) 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


The  advantage  of  our  numerical  approach  is  that  one  can 
calculate  and  tabulate  the  functions  Rjyp,  a;x,  and  Imfii'd 
for  a  given  configuration,  determined  by  the  focal  length  f, 
the  shape  of  the  detection  apertures,  the  width  of  the  deform¬ 
able  optics,  and  the  distance  between  the  crystal  and  the 


detectors.  Then,  to  calculate  the  shape  of  the  interference 
pattern  for  a  certain  phase  distribution  on  the  adaptive  optics, 
one  just  needs  to  change  the  coefficient  of  a  linear  combina¬ 
tion  of  the  tabulated  functions.  This  can  be  a  helpful  tool  for 
studying  the  effect  of  specific  aberrations  on  the  temporal 
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FIG.  3.  (Color  online)  Examples  of  the  shape  of  the  polarization 
quantum-interference  dip  with  two  different  spatial  phase  modula¬ 
tions  (in  black  the  unperturbed  dip,  in  red  the  modulated  one).  In 
the  upper  figure  the  effect  of  a  small  amount  of  coma  along  the 
vertical  axis  is  shown.  In  the  lower  figure  a  more  complicated  su¬ 
perposition  of  aberrations  affects  dramatically  the  shape  of  the  dip. 

interference  or  to  engineer  the  shape  of  the  Hong-0 u-M  andel 
(HOM)  dip. 

Some  examples  of  how  the  temporal  quantum- 
interference  dip  is  modulated  by  a  generic  spatial  phase  shift 
are  reported  in  Fig.  3,  for  coma  (upper  plot)  and  a  superpo¬ 
sition  of  several  different  aberrations  (lower  plot).  The  inter¬ 
ference  visibility  clearly  degrades  in  presence  of  wave-front 
aberrations. 


IV.  PARTICULAR  CASES 

In  this  section  we  will  discuss  the  quantum-interference 
pattern,  described  by  Eq.  (23),  for  a  few  simple  cases.  First 
we  will  consider  the  case  when  no  spatial  modulation  is  as¬ 
signed  to  the  photons  and  Eq.  (23)  will  reduce  to  the  results 
already  described  in  the  literature  for  quantum  interferometry 
with  multiparametric  entangled  states  from  type-1 1  down- 
conversion  [6].  Then  we  will  examine  the  effect  of  a  linear 
phase,  describing  its  implications  for  the  compensation  of  the 
spatial  walkoff  between  the  two  photons.  Finally  we  will 
describe  what  happens  in  the  approximation  of  sufficiently 
large  detection  apertures,  introducing  the  effect  of  even-order 
aberration  cancellation. 

A.  No  phase  modulation 

Applying  no  phase  modulation,  our  equations  reduce  to 
the  ones  derived  in  [6],  Particularly  we  find 


and 


fio=P^(0)PB(0) 


H/G(T)=Sinc 


M\kp  /.  2t 


X 


Mkn 


2diD 


7^2 


Mkn 


2diD 


7^2 


(41) 


(42) 


The  shape  of  the  interference  pattern  is  essentially  given  by 
the  product  of  the  triangular  function  by  the  Fourier  trans¬ 
form  of  the  aperture  function,  centered  at  t=0.  For  physi¬ 
cally  relevant  parameters  the  sine  function  is  almost  flat  in 
the  region  where  the  triangular  function  is  not  zero. 

To  get  an  analytic  result  one  may  assume  Gaussian  detec¬ 
tion  apertures  of  radius  fig  centered  along  the  system's  opti¬ 
cal  axis: 


p(x)  =e 

In  this  case  the  solution  is  quite  simple: 


R(r)=R( 


1-A  1 


DL 


(43) 


(44) 


with 


2diD 

kpMRo' 


(45) 


Typically,  sharp  circular  apertures  are  used  in  experi¬ 
ments.  In  this  case,  the  function  P[q]  is  described  in  terms  of 
the  Bessel  function  J  i(x).  For  a  circular  aperture  of  radius  fi, 


P[q]  = 


ii(2fi|q|) 

fi|q| 


(46) 


However  the  Gaussian  approximation  is  still  a  good  one  if 
the  width  Rq  of  the  Gaussian  is  taken  to  roughly  fit  the 
Bessel  function  (of  width  R):  in  our  case  we  take  fig 

=fi/(2V2). 

Therefore  Eq.  (44)  is  still  approximately  valid  in  the  case 
of  sharp  circular  apertures,  just  taking 


4V2diD 

kpMRg  ' 


(47) 


Mathematically,  in  Eq.  (44)  the  interference  pattern  is 
given  by  the  multiplication  of  a  triangular  function  centered 
at  t=DL/2  by  a  Gaussian  function  centered  at  t=0.  The 
width  of  the  Gaussian  function  ti  decreases  with  increasing 
radius  of  the  aperture  Rb.  Therefore,  in  the  small-aperture 
approximation,  the  width  of  the  Gaussian  is  so  large  that  it  is 
approximately  constant  between  t=0  and  t=DL/2,  giving 
the  typical  triangular  dip  found  in  quantum-interference  ex¬ 
periments.  On  the  other  hand,  increasing  the  aperture  size, 
the  width  of  the  Gaussian  function  decreases,  reducing  the 
dip  visibility  (see  Fig.  4).  Physically,  this  can  be  explained 
by  the  fact  that  by  increasing  the  aperture  size  we  let  more 
wave  vectors  into  the  system,  and  so  the  spatial  walkoff  In 
type-ll  interferometry  introduces  distinguishability,  reducing 


062304-6 


MULTIPARAMETER  ENTANGLED-STATE  ENGINEERING... 


PHYSICAL  REVIEW  A  79,  062304  (2009) 


FIG.  4.  On  the  right  side  we  can  see  the  interference  patterns 
with  three  different  detector  aperture  sizes:  the  corresponding  aper¬ 
ture  functions  are  shown  on  the  left  side. 

the  interference  visibility.  Enlarging  the  aperture  sizes  is  of¬ 
ten  useful  in  practice,  for  example,  to  get  a  higher  photon 
flux.  Moreover,  since  in  the  SPDC  process  different  fre¬ 
quency  bands  are  emitted  at  different  angles,  it  may  be  nec¬ 
essary  to  open  the  detection  aperture  in  applications  where  a 
broader  bandwidth  is  needed.  This  is  clearly  a  problem  when 
using  type-1 1  phase  matching  in  birefringent  crystals,  since 
the  visibility  of  temporal  and  polarization  interference  gets 
drastically  reduced. 

B.  Linear  phase  shift 


FIG.  5.  Quantum-interference  pattern  for  different  detector  ap¬ 
erture  sizes,  introducing  a  linear  modulation  of  the  deformable  mir¬ 
ror,  in  order  to  restore  the  indistinguishability  between  the  photons, 
decreased  by  the  spatial  walkoff  in  the  generation  process. 


In  the  case  of  a  reflective  system,  in  which  the  phase 
modulation  is  implemented  by  means  of  a  deformable  mirror 
(Fig.  5),  tilted  by  an  angle  0, 

<p(x)=2ICotan  0y=Siyy.  (53) 

Therefore,  the  amount  of  tilt  necessary  to  restore  high  vis¬ 
ibility  is 


Suppose  now  we  introduce  a  linear  phase  function  with 
the  spatial  light  modulator,  along  the  direction  Si, 


tan  0  = 


lU 

4f 


(54) 


^(x)=Si-x, 


(48) 


we  get 

Wf/iir)  =Sinc 


2diD  DL 


'  Mk.  . 

'  Mk,  . 

X 

Pb 

In  the  case  of  a  1.5  mm  crystal,  with  M  =0.0723  (pump  at 
405  nm,  SPDC  at  810  nm)  and  lenses  with  focal  length  of  20 
cm  in  the  4f  system,  we  get 

6»  =  0.14  mrad.  (55) 


C .  L  arge-aperture  approximation 


(49) 

If  we  compare  Eq.  (49)  with  Eq.  (44)  we  can  see  that  the 
structure  is  the  same.  We  again  have  a  triangular  function 
centered  at  t=DL/2,  along  with  two  aperture  functions.  But 
this  time,  instead  of  being  centered  at  t=0,  the  aperture  func¬ 
tions  can  be  shifted  at  will  along  the  r  axis.  Suppose  we  now 
apply  a  tilt  along  they  axis  (Si;(=0).  The  modulation  function 
is  then  shifted  to 


'^center  ~  j^/j  ^1/' 

To  get  the  highest  possible  visibility,  the  center  of  the  modu¬ 
lation  function  must  be  matched  to  the  center  of  the  triangu¬ 
lar  dip. 


_DI^ 

“^center  ~  2  ' 


(51) 


SO  that 


If  the  detection  apertures  are  large  enough  for  the  P ,  func¬ 
tion  to  be  successfully  approximated  by  a  delta  function,  we 
get 


W„(r)=  dqG^  -q  Gi  --q  (56) 


Suppose  that  the  spatial  modulator  is  a  circular  aperture 
with  radius  r,  with  unit  transmission  and  phase  modulation 
described  by  the  function  (p(x), 


Gi(x) 


if  |x|  >  r 
if  |x|  <  r. 


(57) 


In  this  case  the  function  <p(x)  can  be  expanded  on  a  set  of 
polynomials  which  are  orthogonal  on  the  unit  circle,  such  as 
the  Zernicke  polynomials: 


V>(q)  =  E  E  ?>nm8n(p)COS(md), 

n  m 
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FIG.  6.  Plot  of  the  values  of  for  /  ,X=-100, ...  ,+100. 
The  radius  of  the  detection  apertures  is  R  =5  mm,  the  distance  be¬ 
tween  the  exit  plane  of  the  modulation  section  and  the  detection 
apertures  is  cf-l=l  m,  and  the  size  of  the  modulation  pixels  is  d 
=25  jubm.  Clearly  only  the  diagonal  elements  are  nonzero,  i.e.,  the 
ones  for  which  X=-l.  In  this  situation  the  effect  of  even-order  ab¬ 
erration  cancellation  is  present. 


m  =  -n,-n  +2,-n  +4, ...  ,n,  (58) 


d  =  0.5  mm 


where  q=(p  cos  0,psin  e).  To  calculate  «p(-q)  we  note  that 
-q=[p  cos(0+77-),p  slnC^+TT-)],  so 

<p(-q)=EEC(p)cos[m(0  +  77)].  (59) 

n  m 

If  m  Is  even  then  cos[m(6'+77)]=cos(m6');  otherwise  If  m  Is 
odd  cos[m(6'+77-)]=  -cosmd).  Therefore 

<p(q)-.p(-q)=2E  2  <p,X(p)cosime).  (60) 

n  m  ocJd 

So,  only  the  Zernicke  polynomials  with  m  odd  contribute 
to  the  shape  of  the  interference  pattern.  This  effect  is  the 
spatial  counterpart  of  the  dispersion  cancellation  effect,  in 
which  only  the  odd-order  terms  In  the  Taylor  expansion  of 
the  spectral  phase  survive.  The  experimental  demonstration 
of  this  effect  was  recently  reported  in  [11]. 

An  interesting  question  is  how  large  the  detection  aper¬ 
tures  should  effectively  be,  in  order  to  obtain  the  even-order 
aberration  cancellation  effect.  According  to  the  numerical  ap¬ 
proach  proposed  in  Sec.  IV,  the  even-order  aberration  can¬ 
cellation  effect  manifests  itself  in  the  limit  where  P(x) 
=  ^(x),  so  that 

a,,,^^(/+\).  (61) 

In  Fig.  6,  a  plot  of  the  value  for  ai  ^  is  shown  for  typical 
values  of  the  relevant  experimental  parameters  (detection  ap¬ 
erture  radius  R  =5  mm,  detection  distance  di=l  m,  and  size 
d=0.1  mm  of  each  pixel  in  the  Fourier  plane  of  the  adaptive 
optical  system).  Clearly,  only  the  diagonal  elements  (the  ones 
for  which  /  =-k)  are  significant,  suggesting  that  the  effect  of 
even-order  aberration  cancellation  may  be  observable  for 
most  typical  experimental  parameters. 

To  get  an  idea  of  what  happens  for  different  experimental 
conditions,  we  can  compute  the  ratio  between  the  intensities 


FIG.  7.  (Color  online)  Plot  of  the  ratio  po  between  the  intensities 
of  the  nondiagonal  coefficient  ^oi  arid  the  diagonal  coefficient  aoo 
as  a  function  of  the  radius  of  Gaussian  detection  apertures,  R,  for 
different  values  of  the  distance  between  the  exit  plane  of  the  modu¬ 
lation  section  and  the  detection  apertures  (di=l,10,100  cm).  On 
the  upper  plot  the  size  of  the  modulation  pixel  d  is  d=0.5  mm, 
while  in  the  second  case  it  is  d=0.25  mm.  Clearly,  for  experimen¬ 
tally  interesting  cases,  the  off-diagonal  coefficient  is  at  least  3  or¬ 
ders  of  magnitude  smaller  than  the  diagonal  one,  leading  to  even- 
order  aberration  cancellation. 

of  the  nondiagonal  coefficient  aoi  ^nd  the  diagonal  coeffi¬ 
cient  aoo: 


The  lower  the  value  for  po  is,  the  less  significant  the  coeffi¬ 
cients  for  \i=-l  are:  the  even-order  aberration  cancellation 
effect  will  therefore  manifest  itself  more  clearly. 

Values  for  po  are  shown  in  Fig.  7  for  two  different  cases. 
In  both  pictures,  the  value  of  po  is  shown  as  a  function  of  the 
Gaussian  detection  aperture  radius  R,  for  three  different  val¬ 
ues  of  the  distance  between  the  plane  113  and  the  detection 
lenses  di.  In  the  upper  panel,  the  size  of  each  small  square  in 
which  the  spatial  phase  is  assumed  to  be  constant  is  d 
=0.5  mm,  while  in  the  second  case  it  is  d=0.25  mm.  In 
both  cases  po  is  significantly  smaller  than  1,  and  it  becomes 
smaller  and  smaller,  increasing  the  value  of  the  detection 
aperture  radius.  However,  po  is  smaller  for  larger  values  of  d, 
implying  that  the  spatial  variability  of  the  modulation  phase 
plays  a  role  in  the  degree  of  even-order  cancellation  of  the 
modulation  itself. 

It  turns  out  that  for  the  aberration  cancellation  effect  to 
appear,  it  is  in  fact  only  necessary  for  one  aperture  to  be 
large  and  for  one  detector  to  be  integrated  over.  This  is  suf¬ 
ficient  to  produce  the  transverse-momentum  delta  functions 
that  lead  to  even-order  cancellation.  To  demonstrate  this,  we 
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can,  for  example,  consider  the  case  where  the  aperture  at  6  is 
large,  and  the  detector  at  B  is  integrated  over,  while  the  ap¬ 
erture  at  A  is  taken  to  be  finite,  with  detector  A  treated  as 
pointlike.  The  location  of  the  pointlike  detector  will  hence¬ 
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forth  simply  be  denoted  asx,  and  we  continue  to  work  within 
the  quasimonochromatic  approximation.  If  we  integrate  only 
over  Xs,  leaving  x  unintegrated,  then  it  is  straightforward  to 
show  that  the  analogs  of  Eqs.  (A3)  and  (A 4)  are 


l/l/^°^(x,q,q', I')  =e  q 


(63) 


l/V(x,q,q',v)  =e  ()(  q  -i- 


aox\ 
cfn  / 


cfn/ 


(64) 


We  have  defined  Q  and  Q'  to  be  the  Fourier  transforms, 
respectively,  of  and  p^.  Pg  is,  as  before,  the  Fourier  trans¬ 
form  of  IpbP-  We  now  let  aperture  B  become  large,  so  that 
the  function  Pg  goes  over  to  a  delta  function.  For  Gi(x) 
=e'<AM  and  G2(x)=l,  we  can  substitute  these  results  into  the 
coincidence  rate  (which  will  now  be  a  function  of  both  rand 
the  position  x  of  detector  A),  and  carry  out  the  q'  and  v 
integrals.  For  the  modulation  term,  we  find 

RM(x,T)=R(x,r)-Ro(x)  (65) 


=j  T/D)e2-qQ(2/q^/kp)[(2T/D)+L] 

'2qh 


xSinc 


kn 


+  0  -q  + 


cfn  / 


cfn 


(66) 


FI  ere,  we  have  used  the  fact  that  the  Fourier  transform  of 
p^(x)  equals  the  complex  conjugate  of  the  Fourier  transform 

of  P/\(-x),  in  order  to  write  Q'  in  terms  of  Q.  We  see  from 
the  presence  of  the  factor  that  even-order  aber¬ 

ration  cancellation  occurs  even  though  one  aperture  is  finite 
and  the  corresponding  detector  is  pointlike.  This  point  may 
be  of  importance  in  future  attempts  to  produce  aberration- 
canceled  imaging. 


V.  CONCLUSIONS 

Summarizing,  in  this  paper  we  have  carried  out  a  theoret¬ 
ical  study  of  the  relation  between  the  wave-front  modulation 
of  the  entangled  SPDC  photons  and  the  shape  of  the  result¬ 
ing  temporal  quantum-interference  pattern.  Due  to  the  mul- 
tiparametric  nature  of  the  generated  entangled  states,  the 
modulation  on  the  spatial  degree  of  freedom  can  affect  the 
shape  of  the  polarization-temporal  interference  pattern  in  the 
coincidence  rate.  Our  aim  is  twofold:  from  one  side  we  want 


to  study  the  effect  of  wave-front  aberration  on  quantum  in¬ 
terferometry,  and  from  the  other  we  want  to  discuss  a  way  to 
engineer  multiparametrically  entangled  states. 

We  have  introduced  a  theoretical  model  for  calculation  of 
the  shape  of  the  polarization-temporal  interference  pattern 
given  a  certain  general  phase  modulation  in  the  crystal  far 
field,  assuming  as  a  free  parameter  the  shape  and  the  dimen¬ 
sion  of  the  collection  apertures.  Using  a  numerical  method  to 
study  the  resulting  equation  has  shown  that  for  typical  ex¬ 
perimental  cases  the  hypothesis  of  large  apertures  can  be 
assumed  to  be  valid.  In  such  an  approximation,  only  the  odd 
part  of  the  assigned  phase  modulation  affects  the  shape  of  the 
interference  pattern.  This  effect  has  recently  been  demon¬ 
strated  experimentally  [11], 

M  oreover,  it  is  often  useful  in  experiments  to  enlarge  the 
collection  aperture  in  order  to  collect  a  higher  photon  flux 
and  larger  optical  bandwidth.  But  when  working  with  type-1 1 
birefringently  phase-matched  down-conversion,  spatial 
walkoff  between  the  emitted  photons  introduces  distinguish- 
ability  between  the  two  possible  events  that  can  lead  to  co¬ 
incidence  detection,  reducing  the  visibility  of  quantum  inter¬ 
ference.  Such  walkoff  can  be  compensated  for  with  a  linear 
phase  shift  in  the  vertical  direction,  restoring  high  visibility. 
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APPENDIX  A:  SKETCH  OF  DERIVATION  OF  EQ.  (23) 


In  this  appendix  we  sketch  the  major  steps  for  the  derivation  of  Eq.  (23).  Substituting  Eq.  (21)  into  Eq.  (16)  and  the  result 
into  Eq.  (22),  one  finds  the  following  expressions  for  Rq  and  I/I/g(t): 

Ro=J  (Al) 

=  cfqQlq'Qlj4)*(q,v)<D(q',-v)G;(^q)Gi(^q')G2(-^q)G2(-^q')l/l/(q,q',v),  (A2) 

where 

l/l/<‘’)(q,q',v)  =  j  dxAdxBH*(XA,q,i')H*(XB,-q,-v)H(XA,q',i')H(XB,-q',-p) 

+  H*(x,,-q,-v)H*(Xs,q,v)H(x,,-q',-v)H(Xs,q',,^)  (A3) 

and 


1/1/ (q 


dxA dxgH  *(X;, , q,  i')H  *(Xb ,  -  q,  -  v)H  (X;, ,  -  q ' ,  v)H  (Xg , q' ,  - 


+  H*(x4,-q,-v)H*(XB,q,i')H  (XA,q',-v)H(XB,-q',v). 
The  angular  and  spectral  emission  function  <I>(q,  v)  is  given  by 


f(q 


■-)=/' 


f 

'z  r 

dzn 

-iA(q,p)z 


Performing  the  integrals  over  the  spatial  coordinates  dx^  and  Xg,  one  gets 

l/l/(^q,q^I^)  =e''(2WI<iP-l‘i'lV/i[(q  -q')]P6[-(q  -q')]  +P/i[-(q  -qOJPsLCq  -q')]} 


(A  4) 

(A  5) 


(A  6) 


and 


l/l/(q,q',  v)  =e'(WI‘'l'-l‘''lV~A[(q  +q')]P6[-(q  +q')]  +P/i[-(q  +q')]PB[(q  +q')]}.  (A  7) 


Finally,  use  of  the  integral  representation  for  the  sine  func¬ 
tion  [Eq.  (A 5)]  allows  the  v  integration  to  be  carried  out,  but 
at  the  expense  of  introducing  two  integrations  over  a  pair  of 
new  parameters  (say,  z  and  z').  Note  the  following  relation, 
which  can  easily  be  verified  by  sketching  the  functions  on 
the  left-hand  side: 


APPENDIX  B:  JUSTIFICATION  OF  EQ.  (32) 

Suppose  we  have  a  set  A,  which  can  be  partitioned  into  a 
collection  of  disjoint  subsets  A*,,  with  k=l,2,...: 

UA|^=A,  Ai(r\Ai  =  (/>  if  k  =A  I .  (Bl) 

k 


n[x]n[x  -  a] 


1  if  -l<a<0,  -|<X<|-Fa 
1  if0<a<l,  -j-I-q;<X<5 
0  Otherwise. 

(A  8) 


From  this,  it  follows  that 


To  each  set  we  can  associate  a  characteristic  function. 


1,  XeA^ 
0,  X  f  Ak, 


(B2) 


such  that 

2  XkM  =xaM,  XkMxiM  =  ^kiXkix) -  (B 3) 

k 


n[x]n[x  -  a]dx  =  A(a) ,  (A  9) 


where  A(a)  is  the  triangle  function.  These  facts  allow  us  to 
carry  out  the  two  z  integrations  that  arise  from  the  sine  func¬ 
tion,  leading  to  the  result  shown  in  Eq.  (23). 


where  xa  is  the  characteristic  function  for  the  full  set. 


AaW  = 


1,  X  e  A 
0,  X  $  A  . 


(B4) 


The  term  assumes  the  value  of  for  ;^(c(x)=l  and 

the  value  of  1  for  A:it(x)=0  [l-;^|^(x)=l],  so 
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exp 


/'E  AXk(x) 

k 


=  n  e''4«w  =  n  {1[1  -  ;t,(x)]  +  e‘<^^Xkix)}  =  11  [1  +  -  Dxkl  ■ 

k  k  k 


If  we  express  the  first  few  terms  we  get 


(B5) 


n  [1  +  -  1)a:J  =  [1  +  (e'^i  -  1)a:i][1  +  - 1)^:2]  •  •  • 

k 

=  1  +  (e'^i  -  l)xi  +  (e'^2  - +  . . .  +  (e'>i  -  i)(e'>2  - +  (e'>i  -  i)(e'>3  - 
+  •  •  •  +  (e'>i  -  l)(e'>i  -  l)(e'>i  -  +  (e'>i  -  l)(e'>2  -  l)(e'>4  -  Dx.xiXa  +■■■■  (B  6) 

So  that  in  the  end 


exp 


(f>kXk(X) 

k 


=  1  +  E  [(e'>k  -  =  1+2  e^'^^xk  -  2  =  2  e'>^. 


Xk- 


(B7) 


Since  the  square  sets  we  have  used  in  Sec.  IV  satisfy  Eq.  (Bl),  then  the  result  expressed  in  Eq.  (B5)  is  valid  for  our  case. 
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We  describe  a  novel  effect  involving  odd-order  dispersion  cancellation.  We  demonstrate  that  odd-  and 
even-order  dispersion  cancellation  may  be  obtained  in  different  regions  of  a  single  quantum  interferogram 
using  frequency-anticorrelated  entangled  photons  and  a  new  type  of  quantum  interferometer.  This  offers 
new  opportunities  for  quantum  communication  and  metrology  in  dispersive  media. 
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Introduction. — The  even-order  dispersion  cancellation 
effect  based  on  nonclassical  frequency-anticorrelated  en¬ 
tangled  photons  has  been  known  in  quantum  optics  for 
some  time  [1,2].  The  nonlinear  optical  process  of  sponta¬ 
neous  parametric  down  conversion  (SPDC)  traditionally 
provides  a  reliable  source  of  frequency-entangled  photon 
pairs  with  anticorrelated  spectral  components,  as  a  conse¬ 
quence  of  energy  conservation.  If  the  frequency  of  the 
signal  photon  is  then  the  frequency  of  its  twin  idler 
photon  must  be  =  Cl p  —  (o^,  where  Clp  is  the  frequency 
of  the  pump  beam.  A  quantum  interferometer  records  the 
modulation  in  the  rate  of  coincidence  between  pulses  from 
two  photon-counting  detectors  at  the  output  ports  of  a 
beam  splitter  in  response  to  a  temporal  delay  between 
two  spectrally  correlated  photons  entering  its  input  ports 
symmetrically.  This  type  of  quantum  optics  intensity  cor¬ 
relation  measurement,  exhibited  in  the  Hong-Ou-Mandel 
(HOM)  interferometer  [3],  is  manifested  by  an  observed 
dip  in  the  rate  of  coincidences.  In  previous  demonstrations 
of  dispersion  cancellation,  one  photon  of  the  down- 
converted  pair  travels  through  a  dispersive  material  in 
one  arm  of  the  HOM  interferometer  while  its  twin  travels 
only  through  air.  The  final  coincidence  interference  dip  is 
not  broadened  in  this  case,  demonstrating  insensitivity  to 
even-order  dispersion  coefficients  [2,4]. 

Even-order  dispersion  cancellation  has  been  used  in 
quantum-information  processing,  quantum  communica¬ 
tion,  and  in  quantum  optical  metrology.  Eor  example,  it 
enhances  the  precision  of  measuring  photon  tunneling  time 
through  a  potential  barrier  [5]  and  improves  the  accuracy 
of  remote  clock  synchronization  [6].  The  same  effect  pro¬ 
vides  superior  resolution  in  quantum  optical  coherence 
tomography  [7]  by  eliminating  the  broadening  of  the  in¬ 
terference  envelope  resulting  from  group  velocity  disper¬ 
sion.  The  potential  of  quantum  even-order  dispersion 
cancellation  has  recently  stimulated  efforts  to  mimic  this 
effect  by  use  of  classical  nonlinear  optical  analogues  [8- 
10]. 

In  this  Letter,  we  introduce  a  novel  type  of  quantum 
interferometer  that  enables  the  demonstration  of  the  odd- 


order  dispersion  cancellation  as  a  part  of  a  new  dispersion 
management  technique.  In  our  design,  both  even-order  and 
odd-order  dispersion  cancellation  effects  can  be  recorded 
as  parts  of  a  single  quantum  interference  pattern. 

HOM  interferometers  are  commonly  used  to  produce 
either  |'T^)  ~  |2,  0)-|0,  2)  state,  when  the  delay  is  set  to 
balance  the  two  paths,  ensuring  destructive  interference  in 
the  middle  of  the  interference  dip,  or  a  superposition  of 
|1,  1),  |0,  2),  and  |2,  0)  states,  when  the  delay  signifi¬ 
cantly  unbalances  two  paths  and  shifts  coincidences  to  the 
shoulder  of  HOM  interference  pattern.  Mach-Zehnder 
(MZ)  interferometers  fed  by  a  particular  quantum  state 
have  also  been  studied  in  detail  [11]. 

In  the  new  design,  two  interferometers  work  together: 
one  output  port  of  a  HOM  interferometer  provides  input  to 
a  MZ  interferometer.  The  state  of  light  introduced  into  the 
MZ  interferometer  is  continuously  modified  when  the 
delay  in  the  HOM  interferometer  is  scanned.  A  signal 
from  one  of  the  HOM  output  ports  is  fed  into  a  MZ 
interferometer  with  a  dispersive  sample  providing  a  phase 
shift  (f)  in  one  arm,  as  shown  in  Eig.  1.  The  delay  T2  inside 
the  MZ  interferometer  is  kept  at  a  fixed  value.  A  peculiar 
quantum  interference  pattern  is  observed  in  the  rate  of 


FIG.  1  (color  online).  Schematic  diagram  of  the  optical  setup. 
The  SPDC  source  produces  pairs  of  frequency-anticorrelated 
photons  combining  on  a  beam  splitter  in  a  HOM  interferometer 
configuration.  Photons  exiting  one  HOM  port  are  fed  into  a  MZ 
interferometer.  Coincidence  events  are  registered  between  two 
single-photon  detectors  at  the  output  ports  of  the  MZ  interfer¬ 
ometer.  A  dispersive  sample  in  one  arm  of  the  MZ  interferometer 
generates  a  phase  delay  (</>). 
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coincidences  between  two  photon-counting  detectors,  Dl 
and  D2,  at  the  output  ports  of  the  MZ  interferometer  as  a 
function  of  Ti.  The  interference  profile  has  two  distinct 
patterns.  The  central  interference  pattern  depends  only  on 
even-order  dispersion  coefficients,  while  the  peripheral 
pattern  depends  only  on  odd-order  terms.  This  ability  to 
manipulate  and  evaluate  odd-order  and  even-order  disper¬ 
sion  terms  independently  in  a  single  quantum  interferome¬ 
ter  opens  new  perspectives  in  quantum  communication  and 
in  precise  optical  measurement. 

Theoretical  model — Eor  detectors  Di  and  D2,  much 
slower  than  the  temporal  coherence  of  the  down-converted 
photons,  the  coincidence  rate  in  such  intensity  correlation 
measurements  is  [12] 

Rci^h  T2)  =  j dti  j  dt2G^^Hti,  t2),  (1) 

with  12)  the  second-order  correlation  function 

12)'- 

t2)  =  K0|#+Vi)4^Hl2)l^)P-  (2) 

and  are  the  electrical  field  operators  at  the 

surfaces  of  detectors  and  D2,  respectively. 

^  / do}je-‘‘^djbj{coj),  (3) 

where  bj{(Oj)  is  the  mode  operator  at  detector  7,  expressed 

in  terms  of  the  input  field  operators  dj{cOj)  [12].  The 
quantum  state  of  light  emitted  in  a  frequency-degenerate 
noncollinear  type-I  phase-matching  SPDC  process  with  a 
monochromatic  pump  Tip  is 

1^)  oc  j da)fico)dl(D,o  +  -  «)|0),  (4) 

where  f{(o)  is  a  photon  wave  packet  spectral  function 
defined  by  the  phase-matching  condition  in  the  nonlinear 
material,  Oq  =  Tlp/2  is  a  central  frequency  of  each  wave 
packet,  =  Oq  +  is  the  signal  photon  frequency,  and 
=  flo  “  ^  is  the  idler  frequency. 

The  phase  shift  (pico)  acquired  by  the  broadband  optical 
wave  packet  as  it  travels  through  a  dispersive  material 
could  be  expanded  in  a  Taylor’s  series  [13]: 

=  Co  +  Clico.i  -  Oo)  +  C2ioj,i  -  Elof 

+  +  ■  ■  ■  ,  (5) 

where  the  linear  term  Ci  represents  the  group  delay  and  the 
second-order  term  C2  is  responsible  for  group  delay  dis¬ 
persion.  In  a  conventional  white-light  interferometer,  Ci  is 
responsible  for  a  temporal  shift  of  the  interference  pattern 
envelope,  C2  causes  its  temporal  broadening,  while  C3 
provides  a  nonsymmetric  deformation  of  the  wave  packet 
envelope.  Higher-order  terms  might  be  included  when  a 
strongly-dispersive  material  is  used  or  in  the  case  of  ex¬ 
tremely  broadband  optical  wave  packets. 


In  the  optical  setup  of  Eig.  1,  the  dispersive  material 
providing  phase  shift  (p{co)  could  be  situated  in  three 
possible  locations.  When  the  sample  is  placed  on  an  arm 
of  the  HOM  interferometer,  it  leads  to  the  well-known 
even-order  dispersion  cancellation  effect  [4].  It  may  be 
shown  that  the  presence  of  a  dispersive  material  between 
the  two  interferometers  does  not  affect  the  coincidence 
interferogram.  We  thus  concentrate  on  the  most  interesting 
case:  we  place  the  dispersive  sample  of  phase  shift  cpico) 
inside  the  MZ  interferometer,  with  delay  T2  set  to  a  fixed 
value,  and  ti  as  the  variable  parameter. 

Eollowing  the  usual  formalism  [12],  one  can  show  that 
the  coincidence  rate  between  the  detectors  is 

Rci^h  Ti)  =  J  d(o[^o  -  T2)  -  T2)] 

X  [/(«)/*(«)  +  f(co)r(-co)e-^‘‘^H  (6) 
where  Oq  is  a  constant, 

and 

T2)  =  + 

Although  not  obvious  from  the  form  of  Eq.  (6), 
is  a  real  function  for  any  spectrum  f{(o),  as 
can  be  seen  by  rewriting  Eq.  (6)  in  manifestly  real  form: 

Tl)  =  j  d(o{\f{o))\^  +  |/(-w)P 

+  [e-2--i/(a;)r(-«)  +  c.c.]} 

X  [‘*>0  -  (9) 

This  fact  ensures  that  the  technique  demonstrated  here 
applies  to  all  types  of  broadband  frequency-anticorrelated 
states  of  light,  including  those  with  nonsymmetric  spectral 
profiles  produced  in  chirped  periodically-polled  nonlinear 
crystals. 

The  final  coincidence  counting  rate  Rc{ti,  T2)  of  Eq.  (6) 
may  also  be  written  as  a  linear  superposition: 

T2)  =  B  +  Ro(Ti)  -  Reven(Th  Tl)  “  RoddiTu  Tl)- 

(10) 

The  first  coefficient  B  incorporates  all  terms  that  are  not 
dependent  on  the  variable  delay  Ti,  providing  a  constant 
after  integration.  It  establishes  a  baseline  level  for  the 
quantum  interferogram.  The  following  terms, 

^o(^i)  =  4  J  dcofico)ri-co)e-^“^^^,  (11) 

^evenC'^l.  ^2)  =  J  d(of((o)f*(-(o) 

.  ^-2/wTi|^^-2/floT2^-/0(flo-w)^-/0(Oo  +  w) 

+  ^2/floT2^/0(Oo-w)^/0(Oo  +  w)j^  (^12) 
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^odd('^i.  Tl)  =  j  d CO f(co)f*(-co) 

.  ^^-2io){Ti  +  T2) ^i(f>iflo-o)) ^-i(f){rLQ  +  co) 

+  ^-2ico{Ti-T2)^-i(f>i^lo-co)^i(f){nQ  +  co)-^ 

are  responsible  for  the  shape  of  the  interference  pattern. 

The  term  /^o('^i)  represents  a  peak  centered  at  =  0 
that  is  simply  a  Fourier  transform  of  the  down-converted 
radiation  spectrum  and  is  insensitive  to  the  dispersion 
associated  with  (pico).  Since  T2)  is  dependent  on 

the  sum  </>(flo  ~  oj)  +  it  is  sensitive  only  to 

even-order  terms  in  the  expansion  Eq.  (5).  This  manifests 
odd-order  dispersion  cancellation  and  generates  a 
dispersion-broadened  function  centered  around  =  0. 
The  last  term  '^2)^  ir^  contrast,  is  sensitive  only  to 

odd-order  dispersion  terms  in  ^(co).  This  term  demon¬ 
strates  the  well-known  even-order  cancellation.  The  coef¬ 
ficients  ^-2iw(ti+T2)  ^-2i(o(ti-T2)  (Jjpg  away 

from  the  center  of  the  interference  pattern  in  opposite 
directions.  Such  decomposition  of  quantum  interference 
terms  makes  it  possible  to  observe  odd-order  and  even- 
order  dispersion  cancellation  effects  in  two  distinct  regions 
of  the  coincidence  interferogram. 

Example. — Our  results  are  illustrated  by  a  numerical 
example  of  quantum  interference  for  a  3 -mm  thick  slab 
of  a  strongly-dispersive  optical  material  ZnSe,  inserted  in 
one  arm  of  the  MZ  interferometer  to  provide  the  phase  shift 
(f){(o).  In  this  experiment,  we  assume  the  use  of  frequency- 
entangled  down-converted  photons  with  a  100-nm  wide 
spectrum.  As  illustrated  in  Fig.  2,  one  can  identify  the 
narrow  peak  Ihe  center,  which  is  insensitive  to 

dispersion,  along  with  the  component  /^even('^n  ^2)^  which 
is  broadened  by  even-order  dispersion  contributions  only. 
This  central  component  of  the  interferogram  illustrates  the 
odd-order  dispersion  cancellation  effect. 

Two  symmetric  side  dips  /^odd('^n  '^2)  appear  shifted  far 
away  from  the  central  peak  by  the  group  velocity  delay  Ci 
acquired  by  entangled  photons  inside  the  dispersive  mate¬ 
rial.  However,  this  shift  can  be  controlled  by  properly 
adjusting  the  value  of  the  fixed  delay  T2.  Such  a  simple 
adjustment  moves  both  dips  back  closer  to  the  center  and 
makes  it  convenient  for  observing  both  dispersion  cancel¬ 
lation  features  in  a  single  scan  of  the  variable  delay  line 
(ti)  inside  the  HOM  interferometer  (see  Fig.  2).  The 
appearance  of  asymmetric  fringes  on  the  side  of  two  dips 
is  a  clear  sign  of  the  third-order  dispersion  [13]. 

Discussion. — This  result  can  also  be  understood  physi¬ 
cally  by  analyzing  all  possible  probability  amplitudes  that 
lead  to  measured  coincidence  events  between  D\  and  D2. 
The  MZ  interferometer  input  is  a  pair  of  spectrally- 
entangled  photons  separated  by  time  delay  ;  if  the  lead¬ 
ing  photon  has  a  high  frequency,  the  lagging  photon  will 
have  a  low  frequency,  and  vice  versa.  We  consider  first  the 
case  when  no  dispersive  element  is  present,  so  that  the  MZ 
interferometer  introduces  only  a  time  delay  T2  between  its 


FIG.  2  (color  online).  The  normalized  coincidence  rate  as  a 
function  of  Tj  when  a  3 -mm  thick  ZnSe  sample  is  placed  in  the 
MZ  interferometer.  The  fixed  delay  T2  =  26  ps  is  used.  The 
insert  illustrates  the  odd-order  dispersion  contribution. 

two  arms.  We  assume  that  T2  is  much  greater  than  the 
photon  wave  packet  width,  r^.  To  explain  the  dependence 
of  the  photon  coincidence  rate  on  ti  ,  as  shown  in  Fig.  2,  we 
consider  three  processes  occurring  at  the  input  ports  of  the 
last  beam  splitter  in  the  MZ  interferometer:  (1)  If  |  ti  |  > 
and  |t2  —  Til  >  T^,  then  the  two  photons  arriving  at  the 
final  beam  splitter  will  be  distinguishable,  so  that  no  quan¬ 
tum  interference  is  exhibited.  (2)  If  IrJ  ~  IT2I,  so  that 
\t2~'Ti\^Tc^  then  quantum  interference  can  occur 
when  the  leading  photon  takes  the  long  path  of  the  MZ 
interferometer  and  the  lagging  photon  takes  the  short  path. 
The  two  arrive  almost  simultaneously  (within  a  time  r^)  at 
the  two  ports  of  the  final  beam  splitter.  Then  the  HOM 
effect  is  exhibited  at  the  beam  splitter,  albeit  with  only  25% 
visibility  because  of  the  presence  of  the  other  possibility 
that  both  photons  arrive  at  a  single  port,  leading  to  a 
background  coincidence  rate  independent  of  Ti.  From  a 
different  perspective,  one  may  regard  this  scenario  as 
similar  to  that  obtained  in  a  Franson  interferometer  [14], 
for  which  photon  pairs  follow  long-long  or  short-short 
paths.  This  scenario  explains  the  components  of  the  coin¬ 
cidence  interferogram  near  =  ±T2,  and  in  this  case  the 
two  spectrally-entangled  photons  entering  separate  ports  of 
the  final  beam  splitter  lead  to  quantum  interference  accom¬ 
panied  by  even-order  dispersion  cancellation.  (3)  Finally, 
when  I  Til  <  T^,  then  one  possibility  is  that  the  photons 
arrive  at  separate  input  ports  of  the  final  beam  splitter. 
Since  these  photons  are  separated  by  a  time  T2  »  r^, 
they  are  distinguishable  and  do  not  contribute  to  quantum 
interference.  The  other  possibility  is  that  the  pair  arrive  at 
the  same  beam  splitter  input  port.  In  this  case,  upon  trans¬ 
mission  or  reflection  at  the  beam  splitter  there  are  two 
alternatives  for  producing  coincidence:  transmission  of 
the  high-frequency  photon  and  reflection  of  the  low- 
frequency  photon,  or  vice  versa.  This  explains  the  compo- 
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nent  of  the  coincidence  interferogram  near  ~  0.  In  this 
scenario,  which  involves  two  spectrally-entangled  photons 
entering  a  single  port  of  a  beam  splitter,  quantum  interfer¬ 
ence  is  accompanied  by  odd-order  dispersion  cancellation. 
We  thus  see  that  the  quantum  interference  effects  exhibited 
in  scenarios  (2)  and  (3)  are  accompanied  by  dispersion 
cancellation — although  in  opposite  manners  in  the  two 
cases. 

In  conclusion,  we  have  demonstrated  a  new  effect  in 
which  even-  and  odd-order  dispersion  cancellations  appear 
in  different  regions  of  a  single  interferogram.  This  is 
achieved  via  frequency-anticorrelated  photons  in  a  new 
quantum  interferometer  formed  by  a  variable  delay  HOM 
interferometer  followed  by  a  single-input,  fixed-delay 
Mach-Zehnder  interferometer.  The  possibility  of  indepen¬ 
dently  evaluating  even-  and  odd-order  dispersion  coeffi¬ 
cients  of  a  medium  has  potential  for  applications  in 
quantum  communication  and  in  quantum  metrology  of 
complex  dispersive  photonics  structures.  In  particular,  the 
ability  to  accurately  characterize  higher-order  dispersion 
coefficients  is  of  great  interest  in  the  study  of  flattened- 
dispersion  optical  fibers  [15,16]  and  in  dispersion  engi¬ 
neering  with  metamaterials  [17].  The  demonstrated  poten¬ 
tial  of  even-order  dispersion  cancellation  has  stimulated 
the  search  for  classical  analogues  [8,9].  We  expect  that  the 
scheme  presented  here  would  also  trigger  the  similar  de¬ 
velopment  of  nonlinear  optical  techniques  mimicking  this 
quantum  effect.  Finally,  note  that  our  apparatus  may  be 
extended  by  adding  a  second  Mach-Zehnder  to  the  unused 
HOM  output  port,  allowing  the  investigation  of  new  four- 
photon  interference  effects. 
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We  investigate  cancellation  of  spatial  aberrations  induced  by  an  object  placed  in  a  quantum  coincidence 
interferometer  with  type-II  parametric  down  conversion  as  a  light  source.  We  analyze  in  detail  the  physical 
mechanism  by  which  the  cancellation  occurs  and  show  that  the  aberration  cancels  only  when  the  object  resides 
in  one  particular  plane  within  the  apparatus.  In  addition,  we  show  that  for  a  special  case  of  the  apparatus  it  is 
possible  to  produce  simultaneous  cancellation  of  both  even-order  and  odd-order  aberrations  in  this  plane. 
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I.  INTRODUCTION  AND  BACKGROUND 
A.  Introduction 

Aberration  or  spatial  dispersion  occurs  when  light  passing 
through  or  reflecting  off  of  an  object  gains  unwanted  phase 
shifts  that  vary  in  the  transverse  spatial  direction  (orthogonal 
to  the  optical  axis).  These  phase  shifts  are  “unwanted”  in  the 
sense  that  they  differ  from  those  obtained  from  Gaussian 
optics  and  cause  distortions  of  the  outgoing  wavefronts. 
Mathematically,  we  can  represent  the  aberrations  by  pure 
imaginary  exponentials  where  x  is  the  transverse  dis¬ 

tance.  Often  (/>(x)  may  be  expanded  into  a  power  series  in  |x| 
and  separated  into  even  and  odd  orders, 


</>(x)  =  +  </>orfrf(x), 

(1) 

<t>eveni'^)  =  S  ayd^Pyid), 

j 

(2) 

(f>even(^)  =  S  Pyi^)  ■ 

(3) 

j 


Here,  r=|x|,  while  P2j{0)  and  P2j+\{0)  are  polynomials  in 
sin  0  and/or  cos  9.  Usually,  the  expansion  is  expressed  in 
terms  of  Seidel  or  Zemike  polynomials  [1-3],  but  for  our 
purposes  the  details  of  the  expansion  are  not  important.  The 
important  point  here  is  simply  that  the  even-order  terms  are 
symmetric  under  reflection,  =  while  the 

odd  terms  are  antisymmetric,  (f>^^^(x)  =  -(f>odd(-^)- 

In  Refs.  [4,5],  a  particular  type  of  interferometric  device 
was  described,  and  it  was  shown  that  if  an  object  was  placed 
in  either  arm  of  this  device,  then  all  even-order  phase  shifts 
introduced  by  the  object  will  cancel  in  a  temporal  correlation 
experiment.  The  effect  is  very  similar  to  the  even-order 
frequency-dispersion  cancellation  first  described  in  Refs. 
[6,7].  As  a  light  source,  the  aberration-cancellation  experi¬ 
ment  used  photon  pairs  produced  via  spontaneous  parametric 
down  conversion  (SPDC).  The  cancellation  effect  depended 
on  the  entanglement  of  the  transverse  spatial  momenta  in  the 
resulting  entangled-photon  pairs. 

In  this  paper  we  re-examine  the  setup  of  Refs.  [4,5]  with 
two  purposes  in  mind.  After  reviewing  the  apparatus  and  the 
even-order  aberration-cancellation  effect  in  the  next  subsec¬ 
tion,  we  first  show  (in  Sec.  II)  that  for  a  special  case  of  the 
apparatus  we  can  in  fact  cancel  all  aberration,  both  even 


order  and  odd  order.  This  cancellation  only  occurs  when  the 
sample  is  placed  in  one  particular  plane  and  opens  up  the 
possibility  of  cancelling  sample-induced  aberration  in  dy¬ 
namic  light  scattering  [8,9],  fluorescence  correlation  spec¬ 
troscopy  [10],  or  other  temporal  correlation-based  experi¬ 
ments.  Our  second  purpose  (carried  out  in  Sec.  Ill)  is  to 
analyze  in  more  detail  the  results  for  the  coincidence  rate  in 
order  to  better  understand  the  physical  mechanisms  involved 
in  aberration  cancellation.  In  Sec.  IV  we  discuss  the  conclu¬ 
sions  that  can  be  drawn  from  these  results. 

Note  that,  because  we  are  motivated  by  the  desire  to  can¬ 
cel  aberrations,  we  will  use  the  phrase  “aberration  cancella¬ 
tion”  for  convenience  throughout  this  paper,  but  in  fact  we 
mean  the  cancellation  of  all  phase  shifts  arising  in  a  given 
plane  not  just  the  subset  that  differ  from  the  predictions  of 
Gaussian  optics.  In  other  words,  “aberration  cancellation” 
here  means  that  only  the  intensity  of  the  light  is  affected  by 
the  object,  not  the  phase.  So,  for  example,  the  placement  in 
the  object  plane  of  an  ideal  lens,  whose  operation  depends  on 
second-order  phase  shifts,  should  have  no  focusing  power  at 
this  point;  it  will  be  as  if  the  lens  is  not  there. 

B.  Even-order  aberration  cancellation 

Consider  the  setup  shown  in  Fig.  1.  In  the  main  part  of  the 
apparatus,  the  two  branches  each  consist  of  a  Fourier  trans¬ 
form  system  containing  lenses  of  focal  length  /  and  a  sample 
providing  a  modulation  Gj{y)  of  the  beam,  where  7  =  1,2 
labels  the  branch  and  y  is  the  transverse  distance  from  the 
optic  axis.  The  Gj  represent  objects  or  samples  whose  prop¬ 
erties  we  wish  to  analyze,  and  the  goal  is  to  cancel  optical 
aberrations  introduced  by  the  samples.  The  case  where  there 
is  a  sample  only  in  branch  1  is  included  by  simply  setting 
G2=1,  but  we  will  keep  the  more  general  two-sample  case; 
we  will  see  later  that  the  extra  generality  pays  off  by  allow¬ 
ing  useful  additional  effects.  A  controllable  time  delay  r  is 
inserted  in  one  arm  of  the  interferometer.  Since  we  will  be 
referring  to  it  often,  we  give  a  name  to  the  plane  containing 
the  samples,  denoting  this  plane  by  H.  The  H  plane  is  simul¬ 
taneously  the  back  focal  plane  of  the  first  lens  and  the  front 
focal  plane  of  the  second.  The  two  lenses  together  form  a  4/ 
Fourier  transform  system.  We  will  examine  in  a  later  section 
what  happens  when  the  sample  is  moved  out  of  the  H  plane. 
Throughout  this  paper,  we  assume  that  the  sample  is  of  neg- 
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n  plane 


EIG.  1.  (Color  online)  Schematic  view  of  aberration- 
cancellation  setup.  (Distances  and  angles  not  necessarily  drawn  in 
correct  proportions.)  The  horizontally  polarized  signal  travels  in  the 
upper  branch  and  experiences  modulation  Gj,  while  the  vertically 
polarized  idler  experiences  modulation  G2  in  the  lower  branch.  G^ 
and  G2  are  both  located  in  the  plane  H,  halfway  between  the  lenses 
of  focal  length  /.  The  beam  splitter  mixes  the  beams  before  they 
reach  the  detectors  Di  and  D2,  which  are  connected  by  a  coinci¬ 
dence  circuit. 


ligible  thickness  compared  to  all  of  the  other  distances  in¬ 
volved  in  the  apparatus.  We  will  refer  to  the  photon  in  the 
upper  branch  (branch  1)  as  the  signal  and  the  photon  in 
branch  2  as  the  idler.  The  polarizing  beam  splitter  sends  the 
horizontally  polarized  photon  into  the  upper  (signal)  branch 
and  the  vertically  polarized  photon  into  the  lower  (idler) 
branch. 

Photons  are  fed  into  the  system  by  a  continuous-wave 
laser  which  pumps  a  nonlinear  crystal,  leading  to  collin- 
ear  type-II  parametric  down  conversion.  The  frequencies  of 
the  two  photons  are  Hq  —  while  the  transverse  momenta 
are  ±q.  For  simplicity,  assume  the  frequency  bandwidth  is 
narrow  compared  to  Hq-  The  two  photons  have  total  wave 

Oq  it  V  Oq 

numbers  which  will  be  approximated  by  k=—  where 
appropriate.  The  down-conversion  spectrum  is  given  by 


q)(q,  v)  =  sine 


LA(q,  p) 
2 


^/[LA(q,F)]/2 


(4) 


Here,  L  is  the  thickness  of  the  nonlinear  crystal  and  for 
type-II  down  conversion  we  have 

2|qP 

A(q,  v)  =  -  vD-\-  M^2  •  q  + - •  (5) 

^pump 

D  is  the  difference  between  the  group  velocities  of  the  ordi¬ 
nary  and  extraordinary  waves  in  the  crystal,  and  M  is  the 
spatial  walk  off  in  the  direction  02  perpendicular  to  the  inter¬ 
ferometer  plane.  The  last  term  in  A  is  due  to  diffraction  as 
the  wave  propagates  through  the  crystal. 

The  parametric  down-conversion  process  may  be  de¬ 
scribed  by  a  Hamiltonian  of  the  form 

H  =  ihxaldj  +  H.c. ,  (6) 

where  d^  and  di  are  annihilation  operators  for  the  signal  and 
idler  photons.  The  constant  x  includes  the  amplitude  of  the 
classical  pump  field.  Applying  the  time  evolution  operator 


e  to  the  vacuum  state,  we  find  that  the  wave  function 
entering  the  apparatus  from  the  crystal  can  be  written  as 

=  (1  -  I  ’7p/2)|0)  +  10^)  +  97^1^4)  +  •  •  • ,  (7) 

where  rj=xt,  and  \'^2n)  represents  a  term  with  n  photons  in 
the  signal  mode  and  n  in  the  idler  mode.  For  parametric 
down  conversion  we  operate  in  the  regime  where  1 97I  ^  1  so 
that  terms  higher  than  l'T'2)  neglected.  In  addition, 

the  vacuum  term  may  be  ignored  since  it  will  not  contribute 
to  coincidence  detection.  Thus,  effectively  our  wave  function 
is  given  by 

|9p)  «  |'p2)  =  J  dqdv^(q,v)dl(q,D,o+  q,flo-  i^)|0). 

(8) 

Note  that  Gi  and  G2  could  be  produced  by  two  separate 
objects  at  two  separate  points  in  space,  in  which  case  we 
would  need  to  use  a  polarizing  beam  splitter  (PBS)  to  sepa¬ 
rate  the  incoming  beams.  Alternatively,  Gi  and  G2  could 
both  be  produced  by  a  single  object  which  acts  differently  on 
the  two  polarization  states,  in  which  case  we  could  dispense 
with  the  PBS. 

In  the  detection  stage,  two  bucket  detectors  Di  and  D2  are 
connected  in  coincidence.  We  add  adjustable  irises  with  ap¬ 
erture  functions  /?i(xi)  and  P2{^7)  in  front  of  the  detectors. 
We  will  end  up  taking  these  apertures  to  be  of  infinite  width 
but  initially  we  leave  them  in  for  reasons  to  be  explained 
below.  A  lens  of  focal  length  is  placed  one  focal  length  in 
front  of  each  detector.  The  distances  from  the  Fourier  plane 
of  the  main  part  of  the  apparatus  to  the  aperture  and  from  the 
aperture  to  the  lens  are  di  and  ^2-  In  order  to  erase  path 
information  for  the  photons  reaching  each  detector,  a  polar¬ 
izer  at  45°  to  the  polarization  directions  of  both  incoming 
beams  is  placed  in  each  path.  The  two  polarizers  are  oriented 
orthogonal  to  each  other. 

The  full  transfer  function  for  each  branch  is  [5] 

77,«(x„,q7,w)  =  Gjy-q}jHap~a,fiiM,  (9) 
where  the  transfer  function  of  the  detection  stage  is 

Ho  (x„,q7,w)  =  e'«‘'i+‘^2+/D)e-(''<^/2/B)[W2//D)-i]x^g-i(<;i/2^)q," 


Va  is  the  Fourier  transform  of  the  aperture  function, 

-  q^.)  =  I  ,  (11) 

with  a; ={1,2}  labeling  the  detector  and  J={s,i}  labeling  the 
signal  or  idler  branch.  In  these  expressions,  k  is  the  longitu¬ 
dinal  wave  number,  k=^{(ol cY - ^  for  |q|^/:. 

The  nonpolarizing  beam  splitter  mixes  the  incident 
beams,  so  each  detector  sees  a  superposition  of  the  signal 
and  idler  beams.  The  positive-frequency  part  of  the  field  en¬ 
tering  detector  a  is  given  by 


053813-2 


SPATIAL-DISPERSION  CANCELLATION  IN  QUANTUM... 


PHYSICAL  REVIEW  A  80,  053813  (2009) 


^  ci(^cl(x)6  (x)) 

+  q,-,  a))di{qi,  o))].  (12) 

Using  this  field,  we  can  compute  the  amplitude  for  coinci¬ 
dence  detection: 

=  I  cfqdv^{q,v) 


XHi2{x2,- q,-  v)  + 

Xif^i(xi,-q,-  v)H,2(x2,q,p)],  (13) 

where  Hj^{x^,qj,[lQ±  v)  have  been  abbreviated  by 
Hj^ix^,qj,±v). 

The  coincidence  rate  as  a  function  of  time  delay  r  is 


R(t) 


=  J  d^Xid^X2dtidt2lA(xi,X2,ti,t2)i^'  (14) 


As  was  shown  in  [11],  R(t)  will  generically  be  of  the  form 


R(t)=Ro 


where  A(v)  is  the  triangular  function: 

1  -  |x|,  |v|  <  1 


A(x)  =  ■ 


0, 


U  >  1 


(15) 


(16) 


The  T-independent  background  term  Rq  and  r-dependent 
modulation  term  W{t)  were  calculated  in  [5]  to  be: 


/q 


/q 


Ro  =  J  cpqd^q'  sinc[MLe2  •  (q  -  q')]Gll  —  IG' 


XP2(-  q  +  q')e“<'"“'®2(q-q')g2iWi/<;p„„p(q2-q'2)^ 


Wir) 


■kl' 


cfqtfq'  sine 


MLe2-(q  +  q')A\l-^ 


XG 


fq 


/qU  (fq’ 


fq’ 


k  I  k  /  \  k  /  \  k 

XPi(q  H-  q'):P2(-  q  -  q')e-('"'")™2(q-q') 


(18) 


Now  let  the  apertures  be  large,  so  that  the  Vj  become  delta 
functions,  reducing  Eqs.  (17)  and  (18)  to 


R(^ 


I 


d^q 


(19) 


W(t)  = 


j2^g-(2;MT/D)e2q(j» 


fA 

k 


(20) 


Suppose  that  Gj{x)  =  tj{x)e^^J^^\  where  tj  is  real  and  the  ef¬ 
fects  of  aberrations  are  contained  in  the  phase  factor  (pj.  Dis¬ 
regarding  the  background  term  for  the  moment,  we  see  from 
the  presence  in  Eq.  (20)  of  the  factors 


•^5  L-;[<^l(/q/U-<^l(-/q/i)] 

k  I 


(21) 


that  even-order  aberration  terms  arising  from  sample  1  can¬ 
cel  from  the  modulation  term.  The  even-order  aberrations 
from  sample  2  cancel  similarly.  This  is  the  even-order  can¬ 
cellation  effect  of  Refs.  [4,5]. 

It  should  be  remarked  that  the  setup  of  Fig.  1  may  be 
simplified  by  removing  the  lenses  immediately  in  front  of  the 
detectors.  We  have  left  both  the  lenses  and  the  apertures  in 
the  setup  because  together  they  lead  to  the  presence  of  the 
Fourier- transformed  aperture  functions  Vj  in  Eqs.  (17)  and 
(18);  the  delta  functions  that  arise  from  the  Vj  when  the 
apertures  become  large  will  serve  as  convenient  bookkeeping 
devices  in  the  following  sections  as  we  trace  various  terms 
back  to  their  origins.  If  we  choose  to  simplify  the  apparatus 
and  remove  the  lenses,  then  Eq.  (10)  will  be  replaced  by 

Ho  (x„,q,-,  w)  = 

a  J 

(22) 


where  d  is  the  total  aperture-to-detector  distance,  with  corre¬ 
sponding  changes  in  Eqs.  (17)  and  (18).  However,  in  the 
large-aperture  limit  this  does  not  affect  the  coincidence  rate, 
which  will  still  be  given  by  expressions  (15),  (19),  and  (20). 


II.  ALL-ORDER  CANCELLATION 
A.  Aberration  cancellation  to  all  orders 

Now,  consider  the  background  term  Rq  in  Eq.  (19).  It 
depends  on  Gi  and  G2  only  through  the  squared  modulus  of 
each.  Thus  any  phase  changes  introduced  by  or  G2  cancel 
completely;  in  particular,  the  background  term  Rq  exhibits 
cancellation  of  aberrations  of  all  orders  not  just  even  orders. 
In  the  current  situation,  this  Rq  is  of  no  importance,  simply 
being  a  constant  and  having  no  effect  on  the  r-dependence  of 
the  correlation.  However,  the  fact  that  all  orders  of  aberration 
can  be  cancelled  in  the  background  term  raises  the  question 
as  to  whether  it  can  be  arranged  for  this  to  happen  in  the 
modulation  term  as  well. 

It  turns  out  that  the  answer  to  this  question  is  positive:  it  is 
possible  to  use  this  apparatus  to  cancel  all  aberrations  in¬ 
duced  by  a  thin  sample,  of  both  even  and  odd  orders.  The 
means  for  doing  so  is  evident  from  examining  Eq.  (20).  Sup- 
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Nonpolarizing 


/i2(y>x) 


f 

iXfiX(f+z)J 


g;i/2{[y2/(/+j)]+(x2//)}g,i/2[l/(/+0]x"2 


X  ^-«^x"-{[y/(/‘+z)]+(x//)} 


(26) 


y,  x',  x",  and  ^  are  the  transverse  distances  at  the  points 
shown  in  Fig.  2.  The  integrals  can  be  carried  out,  giving  us 
the  result  that 


iXf 


FIG.  2.  (Color  online)  Schematic  view  of  apparatus  in  Fig.  1, 
with  Gi  set  equal  to  G2.  (Distances  and  angles  not  necessarily 
drawn  in  correct  proportions.)  Here  and  G2  are  being  produced 
by  a  single  object.  The  signal  and  idler  are  collinear.  It  is  also 
possible  for  G2  and  G2  to  be  produced  by  two  identical  but  spatially 
separate  objects  interacting  with  noncollinear  signal  and  idler. 


h2(y,x)  =  - e  (ik/2f)[{zx^/f)+2x-y]  _  _  ^  1^28) 

iXf 

So  the  impulse  response  for  one  branch  of  the  apparatus 
from  source  to  Fourier  plane  (not  including  the  detection 
stage)  is 


pose  that  Gi(x)  =  G2(x),  as  shown  schematically  in  Fig.  2. 
This  can  happen  in  one  of  two  ways:  either  two  identical 
samples  may  be  placed  in  the  two  arms  or  it  may  be  arranged 
so  that  the  two  beams  both  pass  through  the  same  sample;  in 
either  case  it  is  necessary  for  the  sample  to  act  in  the  same 
manner  on  both  polarization  states.  The  second  possibility 
will  usually  be  of  more  practical  interest,  since  identical 
samples  will  often  not  be  available.  For  Gi  =  G2,  Eqs.  (17) 
and  (18)  give 


Ro  = 


(23) 


W(t)  = 


^2qe-2iMre2-qlD 


h'A^,x)=  hi{^,y)Gj(y)h2{j,x)(fy, 


(29) 


^(ikl2fi)z(^-A) 


(i^ff 


/ 


(30) 


Fourier  transforming  to  find  the  transfer  function  leads  to 

(31) 


{iklf){x-y)^-{ikl2f){zx^) 


(32) 


Setting  Gi(x)  =  t(x)e^^^^\  we  see  that  all  phases  now  cancel 
from  the  r-modulated  term  W.  Thus,  all  aberrations  induced 
by  the  sample,  of  any  order,  will  completely  cancel  from  the 
coincidence  rate. 


y^^{-ikl2z)y^^-{ikzl2f)x^ 


(33) 


Previously,  for  z=0,  this  transfer  function  was  simply  given 

B.  Condition  for  all-ordercancellation  by 


Up  to  this  point,  we  have  assumed  that  the  objects  pro¬ 
viding  the  modulation  were  located  in  the  plane  labeled  IT  in 
Fig.  1.  Now  we  consider  what  happens  if  the  modulation 
objects  (the  samples)  are  moved  out  of  the  IT  plane  by  some 
distance  Consider  a  single  arm  of  the  apparatus,  as 

shown  in  Fig.  3.  We  will  take  the  distance  z  from  11  to  be 
positive  if  the  sample  is  moved  toward  the  source,  and  nega¬ 
tive  if  moved  toward  the  detector.  Now,  the  impulse  response 
functions  for  the  first  and  second  lens  respectively  in  each 
branch  of  the  system  will  be 


^i(^,y) 


f 

iXfiX(f-z)j 


^ikl2{[y^l(f-zmfimgikl2[  l/(/-z)]x '  ^ 


(25) 


Hj{x,q,o))  =  (constants)  GJ  — 


-iqx 


(34) 


Therefore,  for  z  ^  0,  we  must  make  the  replacement  (up  to 
overall  constants)  (Fig.  3) 

(35) 

in  all  previous  results,  and  Eq.  (9)  now  involves  an  integral 
instead  of  a  simple  product.  (For  z  =  0,  the  factors  in  the  last 
set  of  parentheses  become  proportional  to  leading 

back  to  the  previous  results.)  In  particular,  in  Eqs.  (23)  and 
(24),  the  factor  |Gi(^)p  becomes 
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EIG.  3.  (Color  online)  Blown  up  version  of  a  portion  of  one 
branch  from  apparatus  of  Eig.  1  (or  Eig.  2),  with  the  object  moved 
a  distance  z  out  of  the  central  plane,  11. 


/ 


d^ydyGi{y  +  ^x 


g-(iW2z)(y2-y'2) 


(36) 


Clearly,  the  phase  of  Gj  no  longer  cancels  out  of  this  expres¬ 
sion  since  nothing  forces  y  to  equal  y'.  The  arguments  of  the 
two  factors  of  Gi  are  now  unrelated,  so  that  aberration  can¬ 
cellation  no  longer  occurs. 

So  any  cancellation  that  occurs  can  hold  exactly  only  for 
phases  arising  in  the  IT-plane  of  the  Fourier  transform  sys¬ 
tem.  The  cancellation  is  approximate  in  the  vicinity  of  this 
plane.  For  samples  of  finite  thickness,  the  degree  of  approxi¬ 
mate  cancellation  will  diminish  as  the  thickness  increases. 

Defining  6=y-y',  the  exponential  term  in  Eq.  (36)  be¬ 
comes 


^-{ikl2z){26-y-^) 


(37) 


Assuming  that  Gi(y-6+j^x)  is  slowly  varying  in  e  com¬ 
pared  to  the  variation  of  the  exponential,  we  may  obtain  an 
estimate  of  the  distance  z  over  which  the  sample  may  be 
moved  out  of  the  plane  while  still  maintaining  a  high  degree 
of  aberration  cancellation.  The  aberration  cancels  when  e 
=0,  so  we  may  use  the  maximum  size  of  6  as  a  measure  of 
the  degree  of  failure  of  the  aberration  cancellation.  As  z 
^0,  the  rapid  oscillations  of  the  exponential  term  cause  the 
integral  of  Eq.  (36)  to  go  to  zero,  unless  k\2e‘y-^\  also 
goes  to  zero  at  least  as  fast  as  |z|.  So,  we  must  have 


\2e-y-^\ 


z 

k 


~  |zX|. 


(38) 


From  this,  we  have 


(39) 


where  is  the  maximum  value  of  e.  Let  be  the  maximum 
illuminated  radius  of  the  sample.  Then,  by  requiring  that 
1 6^1  ^  r^,  we  have  the  estimate  that 
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A 

This  is  essentially  a  limit  on  how  far  from  stationarity  we 
may  be  and  still  safely  apply  a  stationary-phase  approxima¬ 
tion.  Actually,  we  may  make  this  limit  a  bit  more  precise. 
Since  two-sample  points  y  and  y'  inside  the  Airy  disk  of  the 
lens  cannot  be  distinguished  from  each  other,  we  may  require 
that  \eM\^Rairy^  where 


1.22/X  ,  , 

Rairy= - ^  (41) 

a 

is  the  radius  of  the  Airy  disk.  By  substituting  this  into  Eq. 
(39),  we  can  thus  conclude  that,  at  most,  the  order  of  mag¬ 
nitude  of  |z|  may  be  given  by 


(42) 


Taking  for  example  the  values  m,  a^l  cm,  / 

—  10  cm,  and  X  — 10“^  m,  this  gives  us  an  upper  limit  of 
about  1  mm. 


C.  Comparison  with  dispersion  cancellation 

The  idea  of  aberration  cancellation  via  entangled-photon 
interferometry  arose  in  analogy  to  the  similar  dispersion- 
cancellation  effect  [6,7].  It  is  known  that  even-order  and 
odd-order  dispersion  effects  may  be  separated  so  that  either 
even-order  terms  or  odd-order  terms  may  be  cancelled  [12] 
but  that  it  is  impossible  to  simultaneously  cancel  both  sets  of 
terms  together.  Thus,  it  is  a  surprise  that  in  the  case  of  aber¬ 
rations  such  a  simultaneous  cancellation  should  be  possible. 

The  fact  that  aberration  cancellation  only  occurs  in  a 
single  plane  sheds  some  light  on  the  difference  between  ab¬ 
erration  cancellation  and  dispersion  cancellation.  Aberrations 
are  caused  by  phase  differences  between  different  points  in  a 
plane  transverse  to  the  propagation  direction  of  the  light, 
while  dispersion  comes  about  as  a  result  of  phase  differences 
accumulating  along  the  propagation  direction.  We  have  man¬ 
aged  to  cancel  all  orders  of  aberration  produced  by  a  single 
transverse  plane.  But  since  dispersive  effects  accumulate 
longitudinally,  we  cannot  arrange  their  cancellation  in  all  of 
the  infinite  number  of  transverse  planes  the  photon  travels 
through;  thus,  although  even-order  and  odd-order  dispersion 
may  each  occur  separately,  simultaneous  all-order  dispersion 
cancellation  will  not  occur. 

A  more  physical  explanation  can  be  given  for  the  inability 
in  principle  to  cancel  all  orders  of  dispersion.  Suppose  that 
the  index  of  refraction  is  expanded  about  some  frequency  Wq, 

n{(jo)  =  ^0  +  ni{(i>  -  (Uq)  +  ^2^^  ~  +  *  *  *  .  (43) 

The  phase  and  group  velocities  are 

c 

^  n{ci>)  ’ 
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dk 

^  \d(x) 


=  c 


n{(j())  +  (o 


dn((x)) 

d(jL> 


=  c[nQ  +  2^i(a)  -  (Dq)  +  3^2(^  “  +  *  •  •]  ^ 

(45) 


If  both  the  odd-order  and  even-order  dispersion  coefficients 
vanish  simultaneously  (including  the  zeroth-order  term),  then 
and  ^  both  vanish.  In  consequence,  the  phase  and 
group  velocities  both  diverge.  This  is  in  contradiction  to  spe¬ 
cial  relativity,  which  requires  a  finite  group  velocity.  In  con¬ 
trast,  no  similar  obstacle  exists  to  prevent  the  spatially  dis¬ 
tributed  phase  shift  (/>(x)  from  vanishing,  so  there  is  no 
fundamental  principle  preventing  all-order  aberration  cancel¬ 
lation. 

One  further  point  to  note  is  that  the  dispersive  and  aber- 
rative  cases  considered  here  are  not  entirely  analogous,  in  the 
sense  that  one  is  not  simply  obtained  from  the  other  by  in¬ 
terchanging  time  and  space.  In  the  aberration  case,  the  phase 
is  a  function  of  the  transverse  position  x  in  the  physical 
coordinate  space.  In  contrast,  for  the  dispersive  case  the 
phase  is  due  to  a /r^^^/^^cy-dependent  index  of  refraction; 
i.e.,  the  source  of  the  effect  is  in  the  Fourier  transform 
space,  not  in  the  (temporal)  coordinate  space.  However,  in 
both  cases  the  cancellation  occurs  in  the  Fourier  space.  Thus, 
for  aberration  cancellation  an  optical  Fourier  transform  sys¬ 
tem  is  required  to  move  from  the  coordinate  space  (where  the 
source  of  aberration  is)  to  the  Fourier  space  (where  the  can¬ 
cellation  occurs).  For  the  dispersive  case,  the  source  of  the 
dispersion  already  operates  in  the  Fourier  space  so  it  is  not 
necessary  to  introduce  an  extra  Fourier  transform  via  the 
optical  system. 


III.  PHYSICAL  INTERPRETATION 


We  now  wish  to  develop  a  better  understanding  of  how 
aberration  cancellation  occurs  in  the  polarization-based  coin¬ 
cidence  interferometer  that  we  are  using  to  illustrate  this  ef¬ 
fect.  Let  q  and  q'  be  the  ingoing  and  outgoing  momenta  in 
the  upper  branch  at  the  beam  splitter.  The  ingoing  and  out¬ 
going  momenta  for  the  lower  branch  will  be-  q  and-  q',  as 
in  Figs.  4  and  5  below. 

Note  first  of  all  that  the  coincidence  detection  amplitude 
in  transverse  momentum  space  may  be  written  in  the  form 
A(q)=A^(q)+A^(q),  where  A^  represents  the  amplitude  for 
both  photons  to  be  transmitted  at  the  beam  splitter  and  A^  is 
the  amplitude  for  both  to  be  reflected.  The  counting  rate 
involves  the  integrated  and  squared  amplitude;  if  the  mo¬ 
menta  q  and  q'  were  independent  variables,  we  could  write 
this  as 


j  =  j 


A(q)A'‘(q')d^qd^q' , 


(46) 


which  has  terms  A^(q)A^(q')*+A^(q)A*(q')  involving  inter¬ 
ference  between  reflection  and  transmission  (see  Fig.  4),  as 
well  as  noninterference  terms  A^(q)A^(q')*+A^(q)A*(q') 
(Fig.  5).  However,  q  and  q'  are  not  independent  variables; 
momentum  conservation  and  the  fact  that  the  photons  are 
produced  from  down  conversion  together  force  the  require¬ 


FIG.  4.  (Color  online)  Schematic  representation  of  interference 
terms.  In  the  squared  amplitude  fdqdq'A{q)A^{q'),  the  part  of  the 
amplitude  in  which  both  photons  undergo  reflection  at  the  beam 
splitter  (A^)  interferes  with  the  portion  in  which  both  photons  are 
transmitted  at  the  beam  splitter  (A^).  For  these  terms,  q=-q' ,  due  to 
the  delta  function  that  connects  the  amplitudes. 


ment  q'=±q.  These  constraints  are  explicitly  enforced  in 

the  current  context  by  the  factors  of  Vj  in  Eqs.  (17)  and  (18), 
which  become  delta  functions  in  the  large-aperture  limit.  The 
delta  functions  sew  together  the  amplitudes  A^  and  A^  as 
shown  in  the  figures. 

Suppose  again  that  Gj{x)  =  tj{x)e^^J^^\  Since  we  are  un¬ 
concerned  with  effects  related  to  amplitude  modulation  we 
henceforth  set  tj{x)  =  l.  Examining  Eqs.  (17)  and  (18),  we 
then  note  that  even-order  and  odd-order  aberration  cancella¬ 
tion  arise  from  different  sources.  Even-order  cancellation 
arises  from  the  combination  of  the  following  ingredients: 

(Al)  the  Eourier  transforming  property  of  the  lens  in  the 
focal  plane.  This  converts  the  transverse  momentum  en¬ 
tanglement  into  spatial  entanglement  in  the  H  plane. 

(A2)  The  condition  q=-q'  satisfied  by  the  nonback¬ 
ground  half  of  the  terms  (those  that  comprise  W).  These 


S{q  -  q') 


S{q-q') 


FIG.  5.  (Color  online)  Schematic  representation  of  noninterfer¬ 
ence  terms.  In  the  top  part  of  the  figure  the  transmission  portion  of 
the  amplitude  A^  interacts  only  with  itself,  while  in  the  bottom  part 
the  same  is  true  of  the  reflection  amplitude  A^.  For  these  terms,  q 

=q’ 
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terms  arise  from  the  interference  part  of  the  squared  ampli¬ 
tude,  as  in  Fig.  4. 

(A3)  The  Gy(^)G*(^)  structure  that  arises  from  taking 
the  absolute  square  of  the  amplitude  to  find  counting  rates  in 
quantum  mechanics  (/  =  1 ,2).  Combined  with  the  momentum 
constraint  of  A2,  this  becomes 

In  contrast,  odd-order  cancellation  occurs  when  the  fol¬ 
lowing  combination  of  ingredients  is  present: 

(Bl)  the  Fourier  transforming  action  of  the  lens,  as  in  Al; 

(B2)  for  every  photon  of  transverse  momentum  q  there  is 
a  photon  of-  q  present  due  to  down  conversion;  and 

(B3)  Gi  =  G2  so  that  the  product  Gi(^)G2(-^)  becomes 
=  (Note  that  the  cancellation  is 

taking  place  between  different  terms  of  Eq.  (20)  than  were 
involved  in  the  cancellation  of  A3.) 

In  order  to  have  all-order  cancellation,  there  are  two  pos¬ 
sibilities.  Either  both  of  the  above  sets  of  conditions  may  be 
satisfied  simultaneously,  or  else  a  third  set  of  conditions  may 
be  satisfied: 

(Cl)  same  as  Al  and  Bl; 

(C2)  the  condition  q=q'  must  be  satisfied,  as  in  the  back¬ 
ground  term  this  occurs  in  the  noninterference  terms  of 
Eig.  5;  and 

(C3)  similar  to  A3,  the  Gy(^)G*(^)  structure  arises  from 
the  quantum-mechanical  absolute  squaring  of  the  amplitude. 
But  now,  coupled  with  C2,  we  have  G^(^)G*(^) 
=  giving  cancellation  of  all  orders. 

In  A3  and  C3  the  phase  from  a  single  arm  of  the  interfer¬ 
ometer  cancels  with  itself,  whereas  B3  is  a  cancellation  be¬ 
tween  the  two  different  (but  identical  in  this  case)  arms. 
Cases  A  and  B  both  involve  interference  between  the  ampli¬ 
tudes  and  A^  (shown  schematically  in  Eig.  4),  while  case  C 
comes  from  the  noninterference  terms  of  Eig.  5  and  so  will 
occur  even  if  only  one  of  the  two  amplitudes  A^  and  A^  is 
present. 

IV.  CONCLUSIONS 

To  summarize  the  main  results  of  this  paper,  for  the  ap¬ 
paratus  of  Eig.  1  we  have  found  that: 
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(i)  even-order  aberrations  induced  by  the  samples  Gi  and 
G2  cancel; 

(ii)  if  the  two  beams  overlap  so  that  Gi  =  G2,  then  all 
orders  of  aberration  cancel;  and 

(iii)  these  cancellations  only  occur  if  G^  and  G2  are  con¬ 
fined  to  the  z=0  plane. 

These  results  open  up  the  possibility  of  using  quantum 
interferometry  to  eliminate  the  effects  of  sample-induced  ab¬ 
erration  in  experiments  using  temporal  correlation-based 
methods  such  as  dynamical  light  scattering  or  fluorescence 
correlation  spectroscopy.  Through  the  continued  study  of  ab¬ 
erration  cancellation  and  dispersion  cancellation,  it  is  hoped 
that  a  better  understanding  of  the  effects  of  objects  or  mate¬ 
rials  placed  in  an  optical  system,  and  better  methods  of  con¬ 
trolling  those  effects,  will  gradually  emerge.  The  results  re¬ 
ported  here  are  one  more  step  along  that  path. 

The  effects  described  in  this  paper  make  essential  use  of 
the  spatial  entanglement  (or  equivalently  the  transverse  mo¬ 
mentum  entanglement)  between  the  photons  in  the  down- 
conversion  pair.  In  contrast,  the  frequency  entanglement 
played  no  essential  role.  Similarly,  the  anticorrelation  of  the 
polarizations  was  used  primarily  to  control  the  paths  of  the 
photons  and  then  to  erase  the  path  information;  but  these 
functions  could  be  accomplished  by  other  means.  So  only  the 
spatial  entanglement  was  essential.  On  the  other  hand,  it  is 
the  frequency  entanglement  that  is  essential  for  dispersion 
cancellation.  A  question  for  future  investigation  is  whether 
use  of  the  simultaneous  entanglement  of  frequency,  momen¬ 
tum,  and  polarization  variables  (so-called  hyperentangle¬ 
ment)  may  allow  control  over  further  optical  effects  of  ma¬ 
terials. 
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We  discuss  an  apparatus  that  is  capable  of  producing  correlated-photon  images  that  are  free  of 
object-induced  aberration.  We  show  that  both  quantum-entangled  and  classically  correlated  light 
sources  are  capable  of  producing  the  desired  spatial-aberration  cancellation. 

PACS  numbers:  42.30.Va,42.15.Fr,42.30.Kq 


I.  INTRODUCTION 

Correlated-photon  imaging,  also  known  as  ’’ghost” 
imaging,  in  which  coincidence  measurements  are  used 
to  form  images  via  photons  that  never  interacted  with 
the  object  being  viewed,  has  been  a  topic  of  great  in¬ 
terest  since  its  discovery  using  entangled  photon  pairs 
[l|.  It  has  since  been  found  that  most  aspects  of  ghost 
imaging  can  be  simulated  using  spatially-correlated  clas¬ 
sical  h^ht  J2I,  [j ,  including  thermal  and  speckle  sources 

Separately,  it  has  been  demonstrated  that  the  entan¬ 
gled  photon  pairs  produced  in  spontaneous  parametric 
downconversion  (SPDC)  may  also  be  used  to  cancel  some 
of  the  effects  of  frequency  dispersion  [Elillii  or  spatial 
dispersion  (aberration)  11  [3, 0. 

In  [ll,  it  was  pointed  out  that  it  is  possible  to  con- 
struct  a  device  such  that  if  an  object  is  placed  in  a  par¬ 
ticular  plane  then  all  phase  shifts  induced  by  that  object, 
including  all  orders  of  aberration  induced  by  it,  will  can¬ 
cel.  The  goal  here  is  to  make  use  of  that  observation  for 
imaging.  We  show  that  this  may  be  achieved  by  a  sim¬ 
ple  variation  of  the  traditional  ghost  imaging  apparatus 
of  figure  1.  We  then  show  that,  although  an  entangled 
source  was  required  for  the  temporal  correlation  exper¬ 
iments  discussed  in  [illillil,  only  a  classical  source 
with  transverse  spatial  correlation  is  required  for  imag¬ 
ing. 

The  outline  of  the  paper  is  as  follows.  In  section  II, 
we  review  ghost  imaging  and  aberration  cancellation.  In 
section  III  we  discuss  aberration-cancelled  ghost  imaging 
with  an  entangled  light  source,  followed  by  a  similar  dis¬ 
cussion  with  a  classical  source  in  section  IV.  We  discuss 
an  important  technical  point  about  the  need  for  lenses 
in  front  of  the  detectors  in  section  V,  with  conclusions  in 
section  VI. 


II.  CORRELATED-PHOTON  IMAGING  AND 
ABERRATION  CANCELLATION 

A.  Correlated-photon  imaging 

Correlated-photon  or  ”  ghost”  imaging  [l|  is  done  with 
an  apparatus  like  the  one  depicted  schematically  in  figure 


Correlated 

Photon 

Source 


FIG.  I:  Schematic  depiction  of  correlated-photon  imaging 
setup.  The  photons  in  the  two  arms  have  anticorrelated  trans¬ 
verse  momenta  ±q. 


I.  In  the  original  version,  the  correlated  photon  source 
is  a  nonlinear  crystal  pumped  by  a  laser,  leading 
to  spontaneous  parametric  downconversion.  Entangled 
photon  pairs  with  anticorrelated  momentum  components 
q  and  —  q  transverse  to  the  propagation  direction  travel 
along  the  two  arms  of  the  apparatus.  The  object  to  be 
viewed  is  placed  in  arm  2  (the  upper  branch) ,  followed  by 
a  large  bucket  detector,  D2.  The  detector’s  area  is  inte¬ 
grated  over,  so  D2  can  not  record  any  information  on  the 
position  or  momentum  of  the  photon  that  reached  the 
object;  all  this  detector  is  able  to  tell  us  is  whether  the 
photon  reached  the  detector  unimpeded,  or  whether  its 
passage  was  blocked  by  an  object.  In  the  other  branch  of 
the  apparatus  there  is  no  object,  and  all  of  the  photons 
reach  a  CCD  camera  or  array  of  pointlike  detectors  with¬ 
out  hindrance.  A  coincidence  circuit  is  used  to  record  a 
count  every  time  a  photon  detection  occurs  simultane¬ 
ously  (within  the  coincidence  time  window)  at  each  de¬ 
tector.  By  plotting  the  coincidence  rate  as  a  function  of 
position  xi  in  detector  I,  we  build  up  an  image  of  the 
object.  This  is  true  even  though  photons  that  actually 
encountered  the  object  in  branch  2  left  no  record  of  the 
object’s  position,  and  the  photons  in  branch  I  that  do 
carry  position  information  never  encounter  the  object. 

The  crucial  ingredient  is  the  spatial  correlation  of  the 
downconverted  photon  pair.  The  question  arose  as  to 
whether  the  entanglement  of  the  photons  was  necessary. 
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FIG.  2:  Schematic  of  interferometer  with  even-order 

aberration- cancellation.  Large  bucket  detectors  Di  and  D2 
are  integrated  over. 


or  if  a  classical  source  with  anticorrelated  transverse  mo¬ 
menta  could  mimic  the  effect.  It  was  found  that 

this  was  indeed  possible.  The  correlated  light  source  in 
this  case  consists  of  a  beam  steering  modulator  (a  rotat¬ 
ing  mirror,  for  example)  which  directs  a  classical  light 
beam  through  a  range  of  q  vectors,  illuminating  differ¬ 
ent  spots  on  the  object.  The  beamsplitter  then  turns 
the  single  beam  of  transverse  momentum  q  into  a  pair  of 
beams  with  momenta  q  and  —  q.  The  results  were  simi¬ 
lar  to  those  with  the  entangled  source,  but  with  half  the 
visibility.  It  was  later  shown  that  thermal  and  speckle 
sources  may  also  lead  to  ghost  imaging  (BiSSSi). 


B.  Aberration  Cancellation  in  Quantum 
Interferometry 

Consider  the  setup  shown  in  figure  2  Each 

branch  contains  a  4f  Fourier  transform  system  with  lenses 
of  focal  length  /  and  a  thin  sample  that  provides  modula¬ 
tion  Gj{y)  of  the  beam,  where  j  =  1,  2  labels  the  branch, 
and  y  is  the  position  in  the  plane  transverse  to  the  axis. 
The  goal  is  to  cancel  sample-induced  optical  aberrations 
(position-dependent  phase  shifts  produced  by  the  Gj). 
The  case  where  there  is  a  sample  only  in  one  branch  may 
be  included  by  simply  setting  G  =  1  in  the  other  branch. 
The  plane  containing  the  samples,  which  we  denote  by 
the  letter  11  is  simultaneously  the  back  focal  plane  of  the 
first  lens  and  the  front  focal  plane  of  the  second.  Time 
delay  r  is  inserted  in  one  branch.  The  light  may  be  either 
transmitted  through  each  sample,  or  reflected  off  of  it. 

In  the  detection  stage,  two  large  bucket  detectors  Di 
and  D2^  connected  in  coincidence,  record  the  arrival  of 
photons,  but  not  their  positions.  Apertures  described 
by  aperture  functions  pi(xi)  and  P2(x2)  are  followed  by 
crossed  polarizers  at  45°  to  each  beam’s  polarization,  be¬ 
fore  arriving  at  the  detectors. 

A  continuous  wave  laser  pumps  a  nonlinear  crys¬ 
tal,  leading  to  type  II  parametric  downconversion.  The 
frequencies  of  the  two  photons  are  VLq  with  trans¬ 


verse  momenta  ±q.  For  simplicity,  assume  the  frequency 
bandwidth  is  narrow  compared  to  Qq-  The  two  photons 
have  total  wavenumbers  .  The  downconver- 

c  c 

sion  spectrum  is 


4>(q,  v)  =  sine 


LA(q,  v) 


■  -L  A(q,ty) 

e  2 


where  L  is  the  thickness  of  the  crystal,  and 

2|qP 


A(q,  v)  =  -vD  +  Me2  *  q  + 


(1) 

(2) 


D  is  the  difference  between  the  group  velocities  of  the 
ordinary  and  extraordinary  waves  in  the  crystal,  and  M 
is  the  spatial  walk-off  in  the  direction  62  perpendicular 
to  the  interferometer  plane.  The  last  term  in  A  is  due  to 
diffraction  as  the  wave  propagates  through  the  crystal. 
Ignoring  the  vacuum  term  and  terms  of  higher  photon 
number,  the  wavefunction  entering  the  apparatus  is  ap¬ 
proximately  the  biphoton  wavefunction  of  the  downcon- 
verted  pair,  given  by 


1^')  =  J  (i^qdu  i>(q,j/)a+(q,no  +  J^)a|(-q,Oo-J/)|0), 

(3) 

where  and  hi  are  annihilation  operators  for  the  signal 
and  idler  photons.  For  collinear  pairs,  horizontally  polar¬ 
ized  photons  are  directed  into  the  upper  branch  and  ver¬ 
tically  polarized  photons  into  the  lower  branch  by  means 
of  a  polarizing  beamsplitter.  Alternatively,  noncollinear 
pairs  could  be  used  with  polarizers  selecting  horizontal 
(H)  polarization  in  the  upper  branch  and  vertical  (V)  in 
the  lower  one.  In  either  case,  we  will  refer  to  the  H  pho¬ 
ton  in  the  upper  branch  (branch  2)  as  the  signal  and  the 
V  photon  in  branch  1  as  the  idler. 

The  coincidence  rate  is  of  the  generic  form  [T^ 


R{r)  =  Ro 


1-A  1 


2r 


W{r) 


(4) 


where  A{x)  is  the  triangular  function: 


Hx) 


l-\x\,  |a:|  <  1 

0,  |a;|  >  1 


(5) 


For  large  apertures,  pi(xi)  =  ^2(^2)  ~  1;  so,  as  shown 
in  0  ,  the  background  and  r-modulation  terms  are 


Ro 


W{t) 


(6) 
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d?qe 


^^e2-q 


(7) 


where  k  is  the  longitudinal  wavenumber. 

We  are  concerned  here  with  the  phases  introduced  by 
the  Gj,  not  with  their  effect  on  amplitude,  so  we  will 
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FIG.  3:  Schematic  of  correlated-photon  imaging  setup  with 
aberration- cancellation.  All  orders  of  aberration  cancel. 


assume  that  the  Gj  are  pure  phase  objects,  Gj{x.)  = 
with  aberration  effects  contained  in  the  phase 
factor  (fj.  The  presence  in  equation  (|H|)  of  the  factors 

GJ  (^)  Gi  (— ^)  =  shows  that 

even-order  aberrations  arising  from  sample  1  cancel  from 
the  modulation  term.  The  even-order  aberrations  from 
sample  2  cancel  in  a  similar  manner.  This  is  the  even- 
order  cancellation  effect  demonstrated  in  0  and  0. 
As  pointed  out  in  [lH,  if  Gi  =  G2,  then  both  even- 
and  odd-orders  cancel  simultaneously.  These  cancella¬ 
tions  are  exact  only  for  aberrations  induced  by  passage 
through  thin  objects  in  the  plane  11,  or  by  reffection  from 
a  single  surface  at  the  same  location. 

In  the  term  i?o,  both  even-order  and  odd-order  aber¬ 
rations  occur  even  when  Gi  7^  G2.  For  time  correlation 
experiments,  this  is  an  unimportant  background  term; 
however  it  is  this  term  that  is  the  foundation  of  the  imag¬ 
ing  apparatus  described  below,  since  by  removing  the 
beamsplitter,  the  modulation  term  W{t)  will  be  absent. 

The  physical  mechanism  of  the  various  possible  can¬ 
cellations  are  discussed  in  more  detail  in  [l5[ . 


III.  ABERRATION-CANCELLED  IMAGING 
WITH  ENTANGLED  SOURCE 

Now  we  wish  to  combine  aberration  cancellation  with 
ghost  imaging,  in  the  hybrid  device  of  fig.  3.  In  this 
section,  we  assume  that  the  light  source  is  parametric 
downconversion  from  a  nonlinear  crystal.  This  new  appa¬ 
ratus  differs  from  that  of  fig.  2  in  several  respects.  First, 
we  have  removed  the  time  delay  and  polarization  filters, 
which  are  not  necessary  for  imaging  purposes.  Second, 
the  beam  splitter  before  the  detectors  has  been  removed. 
As  described  in  0,  even-order-only  cancellation  effects 
occur  due  to  interference  between  reffection  and  trans¬ 
mission  amplitudes,  while  non-interference  terms  exhibit 
all-order  interference.  The  beamsplitter  needed  to  pro¬ 
duce  the  interference  effects  desired  in  000  is  not 
needed  for  the  imaging  desired  here.  Without  it  there  will 


be  no  interference,  so  we  expect  to  find  all-order  cancel¬ 
lation.  Finally,  one  bucket  detector  {Di)  is  replaced  by  a 
moveable  pointlike  detector,  or  a  CCD  camera.  Whether 
the  remaining  beamsplitter  is  polarizing  or  not  is  unim¬ 
portant.  After  making  these  changes  and  allowing  for 
an  arbitrary  source  of  correlated  (quantum  or  classical) 
light,  we  arrive  at  an  apparatus  that  looks  very  much  like 
the  ghost  imaging  setup  of  figure  1. 

The  coincidence  rate  at  location  xi  of  Di  is 

R(xi)  =  y*d^X2dtidt2|A(xi,X2,ti,t2)P,  (8) 

where  the  transition  amplitude  is 

^(xi,X2,ti,t2)  =  {0\E['^\xi,ti)E^2^\x2,t2)\^)  ■  (9) 

Taking  the  two  detection  apertures  described  by  pi  and 
P2  to  be  large,  we  can  compute  the  coincidence  rate  to 
be 


R{xi)  =  {[B(xi)  +  B(-xi)]  +  [C'(xi)  +  C^xi)]} 

2 


/ 


where 


Id 


/cxi 


Gi  — xi  G2  — ^xi 


/ 


fo 


B{xi)  =  J  du 


,  (10) 


(11) 


-q.  (12) 


Using  eq.  o,  these  integrals  may  be  evaluated;  they 
turn  out  to  be  xi-independent  constants.  Sweeping  all 
overall  constants  into  a  single  constant  we  find: 


R(xi)  =  TZo 


Gi 


G2 


2 


(13) 


We  see  that  only  the  modulus-squared  of  each  Gj  en¬ 
ters  into  i?(x),  so  that  the  aberration  introduced  by  each 
object  cancels  to  all  orders.  If  G2  =  1,  then  we  have  an 
ordinary  (non-ghost)  image  of  Gi.  On  the  other  hand,  if 
Gi  =  1  then  we  have  an  inverted  ghost  image  of  G2.  In 
either  case,  the  image  is  magnified  by  a  factor  of 


IV.  ABERRATION-CANCELLED  IMAGING 
WITH  CLASSICAL  SOURCE 

We  now  replace  the  downconversion  source  of  the  pre¬ 
vious  section  by  a  classical  source  of  anticorrelated  pho¬ 
tons,  as  in  [ij.  Light  entering  a  beam  splitter  with  trans¬ 
verse  momentum  q  leads  to  outgoing  beams  with  anticor¬ 
related  momenta  q  and  —  q.  If  the  beam  steering  modula¬ 
tor  produces  momentum  spectrum  /(q),  the  input  state 
for  pairs  of  photons  having  the  same  q  before  the  beam¬ 
splitter  will  be  ^  /  d^g'F(q)dJ,(q)dJ,(q)|0),  where  aj,  is  the 
creation  operator  for  pump  photons  and  T(q)  =  /^(q). 
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We  assume  for  simplicity  that  F{q)  is  an  even  function, 
F(q)  =  F(— q).  Denoting  creation  operators  in  the  two 
outgoing  branches  by  a\  and  a^,  the  incoming  photon 
pair  will  produce  a  state  after  the  beamsplitter  given  by: 


1^) 


i  J  (fqF{q)  [al(q) +  4(-q)]  '  [«!  (q)  +  4(-q)]  |0) 

/  (fqF{q)  |^a}(q)a2(— q)  +  additional  termsj  |0),(14) 


where  the  ’’additional  terms”  do  not  contribute  to  coinci¬ 
dence  detection.  The  detection  amplitude  of  eq.  is  then 
proportional  to 


/ 


(fqF{q)e 


zq-(xi-X2) 


ifi(q,Xi)if2(-q,X2),  (15) 


where  Hj  is  the  transfer  function  for  branch  j.  Integrating 
over  1)2,  we  then  have  the  coincidence  rate: 


4^(xi)  = 


k 


fn 


f 


Fl  —XI  ]Gi  ^xi  G2  -^xi 


fn 


f 


fn 


(16) 


This  is  similar  to  the  result  for  the  entangled- source  appa¬ 
ratus,  except  modulated  by  the  factor  F  ^^xi^  which  is 
determined  by  the  details  of  the  beam  steering  modulator. 
Similarly,  for  thermal  or  speckle  sources,  this  factor  will  arises 
from  the  transverse  momentum  spectrum  of  the  source. 


V.  ROLE  OF  THE  DETECTION  LENS 


Consider  now  the  lenses  immediately  before  the  detectors 
in  figure  3.  With  no  such  detection  lens  present,  the  transfer 
function  for  branch  j  would  be 

//,(q,,x,)  =  G,  (17) 


from  which  we  see  that  the  information  from  each  q  value  is 
spread  over  all  x  values.  But  with  the  lens,  eq.  (HU  becomes 


Fj  (qj  7  Xj )  —  Gj 


idicij  2 

e  (5 


(18) 

SO  that  each  q  value  is  localized  at  a  single  point  in  the  de¬ 
tector  plane  via  the  delta  function.  Since  each  q  value  is  also 
matched  to  a  sample  point,  the  localization  in  the  second  case 


defines  a  correspondence  between  points  in  the  sample  plane 
and  points  in  the  detection  plane,  allowing  reconstruction  of 
an  image  by  the  pointlike  detector  Di .  This  can  be  explicitly 
verified  by  computing  the  coincidence  rate  with  or  without 
the  final  lenses,  i.e.  using  either  eq.  dm  or  eq.  dm.  Doing 
this,  we  find  that  without  the  lens  in  branch  1  the  coincidence 
rate  becomes  independent  of  xi,  making  imaging  impossible. 
In  contrast,  removing  the  lens  in  branch  2  has  no  effect.  This 
makes  intuitive  sense:  since  we  integrate  over  X2  anyway,  it 
does  not  matter  whether  the  momentum  information  in  this 
branch  is  localized  or  spread  over  the  entire  detector. 

Thus  we  arrive  at  an  important  technical  point:  the  lens  in 
branch  1  is  essential  for  imaging,  whereas  the  lens  before  the 
bucket  detector  may  be  removed  without  harm. 

The  need  for  a  lens  before  Di  may  be  viewed  as  follows. 
The  4f  system  in  either  branch  transfers  modulation  Gj  from 
the  transverse  coordinate  space  (x)  to  the  Fourier  space  (q), 
which  is  where  the  aberration  cancellation  actually  takes  place 
(see  0).  The  lens  in  front  of  Di  is  then  needed  to  transfer 
the  modulation  back  to  coordinate  space  for  imaging. 


VI.  CONCLUSIONS 


In  conclusion  we  have  proposed  a  means  of  producing  ghost 
imaging  of  a  reflecting  surface  or  thin  translucent  object  with 
cancellation  of  object-induced  aberrations.  The  method  in¬ 
volves  a  relatively  simple  apparatus  and  can  be  done  with 
either  entangled  photon  pairs  or  with  a  classically-correlated 
light  source.  This  aberration-cancellation  effect  has  potential 
for  producing  improved  image  quality  in  a  number  of  appli¬ 
cations  in  biomedical  research  and  other  fields. 
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